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Abstract 



The phenomenon of black hole thermodynamics raises several deep issues 
which any proper theory of quantum gravity must confront: to what extent 
does the inclusion of the back-reaction alter the thermal character of the 
radiation, how can the entropy be understood from a microscopic standpoint, 
what is the ultimate fate of an evaporating black hole, and is the outcome 
reconcilable with unitary time evolution in quantum mechanics? 

In the first part of this thesis, we address the issue of determining what 
the actual emission spectrum from a black hole is, once the gravitational 
field of the emitted quanta is included in a quantum mechanical manner. To 
make the problem tractable, we employ two important approximations: we 
quantize only the s-wave sector of the full theory, and we consider only single 
particle emission. By proceeding in the framework of a Hamiltonian path 
integral description of this system, we are able to integrate out the gravi- 
tational field, thereby obtaining an effective action depending only on the 
matter degrees of freedom. This effective action can then be second quan- 
tized in terms of new, corrected, mode solutions thus enabling the calculation 
of the emission spectrum from modified Bogoliubov coefficients. The results 
are particularly interesting in the case of emission from Reissner-Nordstrom 
black holes, since in the extremal limit our results are dramatically different 
from what a naive, and incorrect, semi-classical calculation would yield. 

The other major topic which we discuss is the dynamics of quantum fields 
on background geometries which undergo quantum tunneling. An example 
of such a system which has important implications for both cosmology and 
quantum gravity in general, is the tunneling of a false vacuum bubble leading 
to the creation of a new universe. To determine what the state of a scalar 
field would be as a result of such a process, we make a WKB approximation of 
the Wheeler-DeWitt equation to obtain an imaginary time evolution equation 
for the scalar field. The state after tunneling is then found by solving this 
equation on a portion of the Euclidean Kruskal manifold, the properties of 
which serve to ensure that at late times thermal radiation emerges. 
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Chapter 1 
Introduction 



Perhaps the most conspicuous deficiency in a physicist's current understand- 
ing of nature is the absence of a usable quantum mechanical theory of gravity. 
Standing in the way has been an utter lack of experimental data to serve as 
a guide toward such a theory. On one hand, the Standard Model of particle 
physics, while not without unresolved issues of its own, provides a coherent 
description of all phenomena which have been observed at existing accelera- 
tors; on the other hand, string theory, while remaining a viable candidate for 
a theory of everything, is presently incapable of making further predictions 
which can be tested by these same accelerators. In the face of such a situa- 
tion there is still reason to believe that progress can be made by extracting 
clues from the structure of the theories that we know today. The remarkable 
quantum mechanical properties of black holes are an important example - 
the phenomenon of black hole thermodynamics is tantalizingly suggestive of 
a more comprehensive theory. 

Before discussing the relevance of black holes to quantum gravity, which 
is the focus of the present work, it will be useful to recall some of the devel- 
opment of the classical theory of black holes. Black holes are hard to see. For 
the experimentalist, this is true for obvious reasons: the defining property 
of black holes as regions of spacetime from which nothing can emerge means 
that their existence can only be inferred, either through their gravitational 
effects on neighboring bodies, or by viewing the characteristic glow of matter 
as it is sucked into a hole. Nevertheless, recent observations have essentially 
removed any doubt as to the existence of black holes in our universe. For 
many years, theorists were similarly unable to see black holes in the equations 



1 



CHAPTER 1. INTRODUCTION 



2 



of general relativity, although in this case the difficulty was due to faulty vi- 
sion rather than a lack of evidence. Although Schwarzschild wrote down the 
metric describing the geometry outside a spherically symmetric body, 



immediately after the discovery of the field equations, an incomplete under- 
standing of the geometry's global structure led to the erroneous conclusion 
that a singularity occurs at the horizon, r = 2M, rendering the solution for 
r < 2M unphysical. Only much later was it realized that the singularity 
was merely a coordinate artifact, and that physical quantities are entirely 
well behaved at the horizon. With this came the understanding that there 
exists a true singularity at r = 0, but it was suspected that deviations from 
spherical symmetry would cause the singularity to be smoothed out in a re- 
alistic collapse process. However, once Penrose |jj proved that singularities 
do occur in classical relativity under very generic conditions, the modern age 
of black hole research was begun. 

Two theorems in classical relativity served to presage the subsequent 
development of black hole thermodynamics. First, it was found that the 
mass and angular momentum of black holes obey the so-called "first law of 
black hole thermodynamics": 



J being the hole's angular momentum, A and f2 the horizon's area and an- 
gular velocity, k is the surface gravity, defined as the force which a person 
at infinity would need to exert in order to keep an object suspended at the 
horizon. For the Schwarzschild hole, k = 1/4M. Second, it was found that 
in any process, the total area of black hole horizons necessarily increases - 
the so-called "second law of black hole thermodynamics" . These results led 
Bekenstein p[ to suggest that a black hole possessed an entropy proportional 
to A, and that the conventional second law of thermodynamics would be up- 
held when this entropy was taken into account. While the formulae suggest 
that something proportional to k plays the role of temperature, this point 
was entirely obscure from the standpoint of classical relativity, since in this 
context a black hole can't emit anything, much less thermal radiation. 



ds 2 = -(1 -2M/r)dt 2 + 



dr 2 



+ r 2 (d9 2 + sin 2 9dtf), 



1 - 2M/r 



dM = —KdA + ttdJ, 
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At this point, quantum mechanics enters the game: Hawking || showed 
that quantum effects cause a black hole to emit exactly thermal radiation 
with a temperature T = k/2tt. This result immediately leads to the identi- 
fication S = A/4. As we shall consider the derivation of this effect in detail 
later, let us simply remark here that it is intimately connected to the exis- 
tence of the horizon as a surface of infinite redshift. This allows particles 
of negative energy to exist inside the horizon, leading to the possibility of 
particle-antiparticle creation, one particle flowing into the hole, one flowing 
out, without violating conservation of energy. Perhaps the most startling 
implication of this result is that, although it only relies on fairly conser- 
vative assumptions about the interaction of quantum fields with gravity, it 
leads to a fundamental conflict with the usual formulation of quantum me- 
chanics. For if the black hole continually emits radiation it will eventually 
disappear altogether, and since thermal radiation is uncorrelated, there will 
not be any trace left of the matter which originally formed the hole. More 
concisely, the process is not described by by unitary evolution since any ini- 
tial state forming the black hole leads to the same thermal density matrix 
after the hole has evaporated. If this is true, then quantum mechanics will 
have to be revised in order to be compatible with this phenomenon of infor- 
mation loss. On the other hand, this conclusion may be premature, as the 
approximation employed in computing the radiance breaks down once the 
hole becomes sufficiently small, so that there is in fact no compelling reason 
why the black hole has to disappear. Various alternatives to the informa- 
tion loss scenario have been the subject of much discussion recently (two 
reviews, written from quite different viewpoints, are j|, ||) but no consensus 
has been established. Much of the recent work has been done in the context 
of two dimensional models which, owing to conformal invariance, allow 
more analytical progress than was previously possible. Unfortunately, the 
key questions remain unanswered. 

The other closely related problem is to obtain a proper understanding of 
the formula S — A/A. It seems quite likely that a resolution of the infor- 
mation loss puzzle is contingent upon understanding in what sense, if any, a 
black hole has the number of states e A ^. One possibility is that a black hole 
has this number of weakly coupled states localized near the horizon, proper 
accounting of which will lead to correlations in the outgoing radiation. Al- 
though the radiation may look approximately thermal, the correlations could 
be sufficient to encode all details of the initial state. Any computation of such 
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correlations must go beyond the free field approximation originally employed 
by Hawking; with this goal in mind, a large part of the present work is 
devoted to obtaining the leading corrections to the free field results. 

In the course of this work we shall be discussing these problems mainly 
in the context of the semiclassical approximation. However, we should point 
out that the validity of this approximation is by no means assured, due 
to the peculiar properties of the horizon. Indeed, it has been argued that 
strong coupling effects [0] and quantum fluctuations || make this description 
unreliable. 

The remainder of this thesis is organized as follows. In Chapter 2 we 
develop the Hamiltonian formulation of gravity, which will play a key role 
in the material that follows. We consider both classical and quantum as- 
pects, paying special attention to the spherically symmetric case. Chapter 
3 provides an overview of quantum field theory in curved space, and a de- 
tailed discussion of black hole radiance. The discussion is phrased in terms 
of a nonstandard coordinate system for the black hole geometry |§, which 
is particularly convenient for doing computations near the horizon. This 
coordinate system also provides insights into the global structure of the ge- 
ometry, as we discuss. Finally, the role of fluctuations in the stress-energy 
is considered, with emphasis on the impact these effects have on the mov- 
ing mirror model. In Chapter 2 we turn to the main focus of this thesis: 
the calculation of self-interaction corrections to the black hole emission spec- 



trum [[!(], plj| . Drawing upon material developed in the previous chapters, 
we show how these corrections can be deduced by calculating the effective 
action for a gravitating thin shell. We then discuss the relation between the 
first and second quantized approaches to this problem, and finally, the diffi- 
culty in obtaining multi-particle correlations by a straightforward extension 
of the single particle calculation. In Chapter 5 we describe some attempts 
at gaining a microscopic understanding of the black hole entropy by count- 
ing fluctuations of the geometry. As we will see, the problem is not how 
to obtain a sufficiently large number of states, but rather how to control 
the divergences which inevitably occur. The topic of Chapter 6 is the effect 
of quantum tunneling upon radiance phenomena PH . The main example 



concerns black hole formation via the tunneling of a false vacuum bubble, a 
process which nicely illustrates the interplay between the Hamiltonian and 
Lagrangian approaches to quantum gravity. We will show how thermal ra- 
diation emerges from the black hole at late times, even though the standard 
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calculation of the radiance is inapplicable. 



Chapter 2 

Hamiltonian Formulation of 
Gravity 

2.1 Classical Theory 

To proceed with the Hamiltonian approach to gravity we perform a 3 + 1 
decomposition of the spacetime manifold, singling out a time coordinate. 
A Hamiltonian is then identified that propagates the three geometry along 
the time direction. The subtlety involved in this procedure is caused by 
the reparameterization invariance of the Einstein-Hilbert action: due to the 
existence of redundant variables, time evolution is not uniquely defined unless 
a gauge is specified. However, there is a well developed formalism for dealing 



with such systems [|13l , as will now be discussed. 

The key to separating out the physical and redundant degrees of freedom 
is to write the metric in 3 + 1 form as |1l4f 



ds 2 = -(N^t) 2 + h ij (dx i + N l dt)(dx j + N j dt) (2.1) 
With this labelling, the Einstein-Hilbert action becomes: 

c G = -r^-V^g k = ^r^hN 1 [ 3 n + k 13 k^ - k\ (2.2) 

l07T l07T 

where 3 7Z is the Ricci scalar associated with hij, K a b is the extrinsic curvature 
of a constant t hypersurface: 

= -^[hij - - N#], (2.3) 
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and K = K%- Throughout, Latin indices are raised and lowered by hij, and 
| denotes covariant differentiation with respect to h^. 

Now, it is seen that no time derivatives of N l or N l appear in the gravi- 
tational action, nor will they when matter is included, so that the canonical 
momenta conjugate to these variables vanish identically: 

dC dC , n JX 

^__0 ; ^ = ^- = 0. (2.4) 

These are constraints on the phase space, h^, on the other hand, have 
nonvanishing canonical momenta n^. These can be computed and used to 
put Cg in canonical form: 

C G = Kij h* j - N'Ht - NiW (2.5) 

where 

Ht = Sirh-^hikhji + tnih jk - h l3 h u )ix l H kl - -^hh 3 IZ, (2.6) 

U l = -27rg. (2.7) 

The constraints Tim = t^n i = must hold at all times, which means that 
their Poisson brackets with the Hamiltonian must vanish. This condition 
yields the secondary constraints 

H t = H i = 0. (2.8) 

All the dynamics of gravity is contained in these constraints. It should be 
stressed that these constraints hold "weakly" — they are to be imposed only 
after all Poisson brackets have been computed. 

2.1.1 Surface Terms 

In deriving the canonical form of the action we ignored all surface terms 
which arose. However, for gravity in asymptotically flat space surface terms 
play an important role and cannot be neglected, so we write 

H G = J (fx^Ht + NiW] + surface terms. (2.9) 
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Regge and Teitelboim |15| showed that Einstein's equations are equivalent 
to Hamilton's equations applied to Hq only for a particular choice of surface 
terms. The point is that if Einstein's equations are written as 



hij = Affair) ; 
then we need the variation of Hq to be 



7T; 



-Bij(h, 7r) 



SH, 



G 



BijSh 13 



(2.10) 



(2.11) 



to ensure consistency with Einstein's equations. To put 8Hq in this form we 
must integrate by parts, because space derivatives of and tt^- appear in 
Hq, and then demand that the resulting surface term cancels with the one we 
have included in ( |2.9| ). In [IS] the surface terms were worked out assuming 
a particular rate of fall-off of h^, n^, N f , and N l at infinity. This condition 
can be relaxed [PJ, and in general the full Hamiltonian is: 



with 



where 



H G = / d^lN'Ht + N{W] + / d Sl [H l ls + H 



l 2s\ 



lO/T 



16tt 



AT* 



VhiWdiN 1 - N l diN*) 



(2.12) 



(2.13) 
(2.14) 



u jk - l jk~ °{j L k)- 

The result of Regge and Teitelboim corresponds to choosing fall-off conditions 
such that H l 2s vanishes at infinity. 

Since TC t and 7i l vanish when the constraints are satisfied, the numerical 
value of the Hamiltonian is given solely by the surface terms. This numerical 
value is the ADM mass of the system. 



2.1.2 Spherical Symmetry 



We now specialize to the case of spherically symmetric geometries |T7, 18), . 
It is possible to proceed substantially further in this case because there are 
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no propagating degrees of freedom - after solving the constraints there will 
only be one free parameter remaining: the ADM mass. The metric is written 

ds 2 = -(N l dt) 2 + L 2 (dr + N r dt) 2 + R 2 (dd 2 + sin# 2 d0 2 ) (2.15) 

and the action is 

— / d^x^gTl = [dtdr \tt r R + ti l L - N l H t - N r H r ] (2.16) 
on J J 



16vr 
with 

JU = ^-™ + (««X ; % — R'ttr — Ln' L (2.17) 
1 2R 2 R \ L J 2L 2 r R L K ' 

where ' represents d/dr. 

As before, the Hamiltonian is 

H = J dr[N t Ht + N r H r ] + M (2.18) 

and the constraints are 

n t = H r = 0. (2.19) 

The constraints can be solved as follows. 7Tr is eliminated by forming the 
linear combination of constraints 

Defining 

.... 7r? R RR! 2 

*M s 2S + 2-W (2 ' 20) 

this constraint is equivalent to M! = 0. By comparing with ( 2.13|) it can be 
shown that Ai(oo) is the ADM mass. The constraints can now be solved for 
the momenta: 



tt l = T]Ry (R'/L) 2 - 1 + 2M/R ; tt r = —/ L (2.21) 



where r] = ±1. 
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Since we have denned the new variable, A4(r), it is natural to ask what the 
momentum conjugate to M.(r) is. From the fundamental Poisson brackets, 

{L(r), vr L (r')} = {R(r), n R (r')} = S(r - r') (2.22) 

it is straightforward to check that 

_ Ly/(R>/L)*-1 + 2M/R 

nM = 1 - 2M/R (2 ' 23) 

satisfies 

{M(r),n M (r')} = 6(r-r'). (2.24) 

ii m has a simple physical interpretation. Consider the geometry in Schwarzschild 
coordinates: 

dr 2 

ds 2 = -(1 - 2M/r s )dt 2 s + - - ' + r 2 s {d6 2 + sin 2 ddf) (2.25) 

/ ^* S 

and then transform to new coordinates t(t s ,r s ), r(t s ,r s ). One finds that in 
the new coordinates, 

R - ^ r > ; 12 - r^MF, (^) 2 - (1 - 2M/r ' ) (2 ' 26) 

which then yields 

*M(r) = ?±. (2.27) 

In other words, tim measures the rate at which Schwarzschild time changes 
along a hypersurface. Pm = / ^M( r )dr is conjugate to the ADM mass: 
{M, P M } = 1- This gives clear expression to the often heard statement that 
the energy generates time translations at infinity. 

Since Pm is related to the behavior at infinity, it is not surprising that it 
is invariant only under small diffeomorphisms: 

5PM = [Jdr[f(r)n t + g(r)n r ],P M ] 

vanishes provided 

/(oo),#(oo) < 

as can be checked explicitly. For most coordinate choices, the inequality 
reduces to /(oo) = g(oo) = 0. 
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2.2 Quantization 

Once the Hamiltonian formulation of gravity has been developed, the tran- 
sition to the quantum theory is quite straightforward. This is true only in 
a formal sense, though, since operator ordering ambiguities will be left un- 
resolved, as will the problem of ultraviolet divergences. Furthermore, the 
interpretation of the resulting theory is not at all clear, as we will briefly 
discuss. More detailed discussion of interpret at ional issues can be found in 

To quantize, we realize the phase space variables hu, 7T{ j clS operators 
satisfying the commutation relations 

[h^x'), Tr mn (x H )] = zc^^V - A (2-28) 
We therefore make the replacements, 

7T J — > -I 



Shij 



The states of the system are then described by functions of the three-geometry, 
i.e. by the wave functionals \P[/iy(a;*)]. Physical states are required to be 
annihilated by the constraints: 

H i ^lh ij ] = 2i(^-) m[h ij \=0 (2.29) 




HMh 



l.h 



2 



m* bhijbhki 167T 



= (2.30) 



where G^i = h^hji + huhjk — h^hki, and we have restored the Planck mass. 
The first set of constraints simply says that physical states should be invari- 
ant under reparameterizing the spacelike hypersurfaces; they are analogous to 
Gauss' law as they are linear in derivatives. The final constraint is quadratic 
in derivatives and is known as the Wheeler-Dewitt equation. In asymp- 
totically flat space, the wavefunction also satisfies a non-trivial Schrodinger 
equation: 

m[hilA = .i?^A. (2.31, 
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It is then assumed that ^[/^-^[/ty] can somehow be interpreted as a prob- 
ability density over three-geometries. 

Even if we knew the correct operator ordering prescription, the full Wheeler- 
Dewitt equation would still be much too difficult to solve. However, it can 
be simplified greatly by employing the WKB approximation, which amounts 
to taking m p —>■ oo. In the WKB approximation the wavefunction is written 
in the form 

■qf[hij] = e im p s[h ^ ] . (2.32) 

Inserting this expression into ( [2.30|) and keeping only terms of order m 2 , we 
find that S satisfies: 

16nG ljM ^^- - ^-hh 3 K = 0. (2.33) 
bhij bhki 167r 

This is the Einstein-Hamilton- Jacobi equation, which means that S[hij] is the 
classical action associated with a solution to Einsteins's equations. Specif- 
ically, given some solution g^x^), S[hij] is found by integrating Sq = 
J dftx^f—g 1Z in the region between some reference hypersurface and the 
hypersurface hij(x % ). Thus the WKB approximation reduces the task of solv- 
ing a functional differential equation to solving a set of partial differential 
equations. 

As we remarked earlier, it is difficult to interpret the wavefunction ^[/ly]. 
In particular, what does it mean to say that |\1/| 2 is the probability to find 
some three-geometry? It may well be that the wavefunction only has meaning 
once matter is included. Then, writing *f?[hij,<f>], we would interpret 



\*[h 



as the relative probability of finding various matter configurations on the 
hypersurface hij. We will return to this intepretation in later sections when 
we couple a scalar field to gravity. 



2.2.1 Spherical Symmetry 

In the spherically symmetric case we have the operators L, R and ttl = 
7r r = —ijp, and the wavefunction, ty[L,R], which satisfies the constraints 



H t V[L,R] = H r ^[L,R] = 0. 



(2.34) 
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In this simplified theory the quadratic constraint is still too difficult to solve 
so we again employ the WKB approximation and write 



*[L,R] 



D iS[L,R] 



(2.35) 



Since the WKB approximation amounts to the replacements ttl 
ss_ 

SS 



J|r, and using (|2.21| ), S is found to satisfy 



6L> n R 



5L 



r,Rj{K/LY 



2M/R 



(2.36) 



8S _ L d 

JR ~ ~R 1 dr~ V'~~ v v_ ~ ' ' J 

These expressions are most easily integrated as follows. Start from some 
arbitrary geometry L, R, and vary L while holding R fixed until (R'/L) 2 — 
1 + 2M/R = 0. The contribution to S from this variation is 

S[L, R]=r]Rj dr dL^J(R'/L) 2 - 1 + 2M/R 



r) J dr 



RLJ(R'/L) 2 -1 + 2M/R 



+ RR' log 



R'/L - J(R'/L) 2 - 1 + 2M/R 



1 - 2M/R\ 



(2.37) 



Now vary L, R, while keeping (R'/L) 2 — 1 + 2M/R = 0, to some set geometry. 
There is no contribution to S from this variation, so (|2.37|) is the complete 
solution. The wavefunction, either oscillates or decays exponentially de- 
pending on the sign of [R! / L) 2 — 1 + 2M/ R. This is also the condition which 
determines whether the hypersurface L, R, can be embedded in the classical 
Schwarzschild geometry. 

The description of the state in terms of the wavefunction ^f[L, R] is, of 
course, highly redundant since, according to the constraints, the wavefunc- 
tion must take the same value for many different L's and i?'s. If we wish to 
eliminate this redundancy we should describe the state in terms of physical 
observables, i.e. those which commute with the constraints. In the present 
case there are only two such observables [19|, ^TJ : the ADM mass, M, and its 
canonical momentum Pm- To satisfy the commutation relation [M, Pm] = i 
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we make the replacement M = id/dPM, so that a state of definite mass is 
written 



It is seen that Pm functions as an intrinsic time variable. 

The emergence of an intrinsic time in this system is quite interesting in 
light of the so-called "problem of time" in quantum gravity. In the asymp- 
totically flat case this problem is relatively benign, since the wavefunction 
depends on the time t, which can be measured by clocks at infinity. How- 
ever, for a closed universe there is no such asymptotic region and no obvious 
variable to play the role of time. This makes the wavefunction of a closed 
universe especially difficult to interpret. In quantum mechanics, we interpret 
the wavefunction ifj(x,t) by saying that at fixed time t the probability to 
find x is \ip(x, t)\ 2 . To carry this procedure over to the quantum theory of 
a closed universe we need to specify a variable which can be fixed in order 
to determine probabilities. In the spherically symmetric, asymptotically flat 
case, Pm, like t, can play this role, but no such variable is known in general. 




(2.38) 




(2.39) 
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This completes our discussion of the quantization of the pure gravity 
theory. In what follows we turn to the effects of incorporating matter into 
the system. 



Chapter 3 



Quantum Field Theory Near 
Black Holes 



3.1 A Useful Coordinate System 

Calculations of physical effects in the presence of black holes are greatly 
simplified by employing appropriate coordinates. In this section, we will 
describe the properties of a little known set of coordinates || which will be 
used repeatedly in the material that follows. 

Schwarzschild found his remarkable exact solution for the geometry out- 
side a star in general relativity quite soon after Einstein derived the field 
equations. Further study of this geometry over the course of several decades 
revealed a series of surprises: the existence and physical relevance of pure 
vacuum "black hole" solutions; the incompleteness of the space-time covered 
by the original Schwarzschild coordinates, and the highly non-trivial global 
structure of its completion; and the dynamic nature of the physics in this 
geometry despite its static mathematical form, revealed perhaps most dra- 
matically by the Hawking radiance ||. Discussions of this material can now 



be found in advanced textbooks |24|], but they are hardly limpid. 

In the course of investigating an improvement to the standard calculation 
of this radiance to take into account its self-gravity, as we detail in later 
sections, we came upon a remarkably simple form for the line element of 
Schwarzschild (and Reissner-Nordstrom) geometry. This line element has an 
interesting history ||22|| , but as far as we know it has never been discussed 
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from a modern point of view. We have found that several of the more subtle 
features of the geometry become especially easy to see when this line element 
is used. 

To motivate the form of the coordinates we reconsider the constraint 
equations of spherically symmetric gravity ( |2.17f) : 

„ «S k l ix r ( RR'Y R' 2 
nt= 2&-^r + \-T) ~ 2L = °' 

H r = R'-kr - Lit' = 0. (3.1) 
The canonical momenta are given by 

N r R' RR (N r LR)' (LR) 

nL = ^~J^ ; 7r * = — iv* W ( ) 

We can arrive at a particular set of coordinates by choosing a gauge and 
solving the constraints. Our choice is simply L = 1, R = r. With this 
choice the equations simplify drastically, and one easily solves to find hl = 
\/2Mr , Ti R = J§ and then N* = ±1, N r = ±y^. Thus for the line 
element we have 



ds 2 = -dt 2 + {dr ± y ^-dt) 2 + r 2 {dd 2 + sin 2 Odcf 2 ) . (3.3) 

M, which appears as an integration constant, of course is to be interpreted 
as the mass of the black hole described by this line element. 

For the Reissner-Nordstrom geometry, the same gauge choice leads to a 
metric of the same form, with the only change that 2M — > 2M — Q 2 jr. 

These line elements are stationary - that is, invariant under translation 
of t, but not static - that is, invariant under reversal of the sign of t. Indeed 
reversal of this sign interchanges the ± in ( |3.3j ), a feature we will interpret 
further below. Another peculiar feature is that each constant time slice dt = 
is simply flat Euclidean space! 

We can obtain a physical interpretation of these coordinates by comparing 



them to those of Lemaitre 23 , in terms of which the line element reads 



ds 2 = -^ 2 + [f( f L ^C)p/3 drl - m) 2l "[\{r L -r)fl\de 2 + ^ 2 



(3.4) 
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As the Lemaitre coordinates are synchronous (g TT = — 1 , g T \ = 0), a class of 
timelike geodesies is given by motion along the time lines (tl, 0, =constant), 
and the proper time along the geodesies is given by the coordinate r. To 
arrive at ( |3.4j) we retain the Lemaitre time coordinate, but now demand that 
the radial coordinate squared be the coefficient multiplying d6 2 + sin 2 8d(p 2 . 
In other words, we write 

t = T ; r = {2Mf^{r L -r)] 2 /\ (3.5) 



A simple calculation then leads from ( |3.4|) to (|3.3| ) with the upper choice of 
sign; the lower choice is obtained by repeating the same steps starting from 
Lemaitre coordinates with the sign of r reversed. 

Finally, let us note that these coordinates are related to Schwarzschild 
coordinates, 



9 /If rlr 2 

ds 2 = -(1 - ^-)dt 2 s + ^ m + r\dd 2 + sin# 2 d0 2 ) (3.6) 



by a change of time slicing, 



t s = t-2V2Mr-2Mlog 



V2M 



(3.7) 



In contrast to the surfaces of constant t s , the constant t surfaces pass smoothly 
through the horizon and extend to the future singularity free of coordinate 
singularities. 

In terms of r and t, the radially ingoing and outgoing null geodesies are 
given by 

ingoing: t + r — 2V2Mr + AM hg[y/r + V2M] = v = constant 
outgoing: t-r- 2v / 2Mr - 4M log [x/r - \FlM\ =u = constant . (3.8) 



3.1.1 Global Structure 

Now let us discuss the global properties of our coordinate system. Perhaps 
the clearest approach to such questions is via consideration of the properties 
of light rays. Taking for definiteness the upper sign in ( |3.3j ), and without any 
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essential loss of generality restricting to the case d9 = d(j) = appropriate to 
the 9 = 7r/2 sections, we find that ds 2 = when 



For the class of light rays governed by the upper sign, we can cover the entire 
range < r < oo as t varies. In particular one meets no obstruction, nor any 
special structure, at the horizon r = 2M. For the class of light rays governed 
by the lower sign there is structure at r = 2M. When r > 2M one has a 
positive slope for ^, and r ranges over 2M < r < oo. When r < 2M one has 
a negative slope for ^, and r ranges over < r < 2M. When r = 2M it does 
not vary with t. From these properties, one infers that our light rays cover 
regions I and II in the Penrose diagram, as displayed in Fig. Let us 

emphasize that the properties of the Penrose diagram can be inferred from 
the properties of the light rays, although we will not belabor that point here. 

If one chooses instead the lower sign in (|3.3|) , and performs a similar 
analysis, one finds that regions I and II' are covered. Patching these together 
with the sectors found previously, one still does not have a complete space- 
time. However our line element is not yet exhausted. For in drawing Figure 
1 we have implicitly assumed that t increases along light rays which point 
up ("towards the future"). Logically, and to maintain symmetry, one should 
consider also the opposite case, that the coordinate t increases towards the 
past. By doing this, one generates coordinate systems covering regions I' 
and II' respectively I' and II', for the upper and lower signs in ( |3.3| ) . Thus 
the complete Penrose diagram is covered with patches each governed by a 
stationary - but not static - metric, and with non-trivial regions of overlap. 

In the Reissner-Nordstrom case the generalization of ( |3.3|) has a coor- 
dinate singularity at r = Q 2 /2M. However this singularity is inside both 
horizons, and does not pose a serious obstruction to a global analysis. One 
obtains the complete Penrose diagram also in this case by iterating construc- 
tions similar to those just sketched. 

The usual Schwarzschild line element appears to be time reversal sym- 
metric, but when the global structure of the space-time it defines is taken 
into account one sees that this appearance is misleading. The fully extended 
light-rays in Figure 1 go from empty space to a singularity as t advances 
(they pass from region I into region II), which is definitely distinguishable 




(3.9) 
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R = 




R = 



Figure 3.1: Penrose diagram for the Schwarzschild geometry. As described 
in the text, r and t in one coordinate patch, for the upper sign of the line 
element, cover regions I and II. As is clear from the diagram, the ingoing 
light rays are captured in their entirety (into the singularity), whereas the 
outgoing light rays cannot be traced back past the horizon. 

from the reverse process. There is a symmetry which relates these to the 
corresponding rays going from region II' to region I', however it involves not 
merely changing the sign of t in the Schwarzschild metric, but rather going to 
a completely disjoint region of the space-time. This actual symmetry of the 
space-time is if anything more obvious in our construction than in the stan- 
dard one. Thus by taking the line-element in region I stationary rather than 
static we have lost some false symmetry while making the true symmetry - 
and its necessary connection with the existence of region I' (constructed, as 
we have seen, by simultaneously reversing the sign of t and interchanging the 
future with the past) - more obvious. 
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3.2 Quantum Field Theory in Curved Space 

In this section we briefly present some of the important results in the theory 
of quantum fields in curved space. Our goal is primarily to establish notation 
and to write down some formulas which will be referred to in later sections. 



Comprehensive reviews of the subject can be found in |24], |25[ . 

Quantum field theory in curved space is a hybrid of quantum and classical 
field theory, which one hopes reliably describes the behavior of matter in 
regions of relatively low space-time curvature. The approach is to couple a 
quantum field to a metric tensor, which is treated as a classical variable. For 
simplicity, we will consider a massless scalar field with action 

S = ~\J d'x^g^d^d^. (3.10) 
The field satisfies the wave equation, 

(V=s<r 0,^ = (3.ii) 

Associated with this wave equation is the conserved inner product 

(0 1; 2 ) = -i J d\^/hrf(j)x{x\ t) <J <j)* 2 (x\ t) (3.12) 

where <pi and 02 are two solutions of the wave equation. Here £ is a Cauchy 
surface, and n M is a future pointing unit vector normal to E. In order to 
second quantize the field we must write down a complete set of "positive 
frequency" solutions, Ui(x^), which satisfy (ui,Uj) = 5ij. This is where the 
ambiguity in the quantization process occurs, since there are many different 
sets of Mj's, which will be shown to lead to inequivalent quantizations. Once 
we have chosen a set, the field operator is written in terms of creation and 
annihilation operators: 

<p = V] [diUi + a\u*} (3.13) 

i 

where a, obey 

[ai,a]} = 5 i:j . (3.14) 

The vacuum state is defined by a.j|0 u ) = 0, and particle states are created by 
applying a^'s to |0 n ). 
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On the other hand, suppose that instead of Ui we had chosen the modes 
Vi in which to expand the field: 

= EM + 6 kl- ( 3 - 15 ) 

i 

Then the vacuum would be defined by bi\O v ) = 0, and particle states would 
be created by applying &t's to |0„). The obvious question which arises is: 
what is the relation between the states defined by Ui and those defined by 
Vi. To answer this, we note that by equating fl3.13| ) and ( |3.15| ), and taking 
inner products, we obtain the "Bogoliubov transformation" : 

&i = E(«>* - £>!) ( 3 - 16 ) 

where 

a ij = (v i ,u j ) ; Pij = -(vi,Uj) (3.17) 
The Bogoliubov coefficients obey the completeness relations, 

k 

J2( a ikPjk ~ PikOijk) = 0. (3.18) 

k 

It can then be shown that the states are related by 

|V>«) = C :exp \^a(a~ l (3)a + a{a^ - l)a ] + ]^a ] (- (3* a' 1 )^} : \if> v ) (3.19) 
where C is a constant. The average number of v particles in the u vacuum is 

(oj&J&iioo = Y,\M- ( 3 - 2 °) 



3.3 Black Hole Radiance 

Now we turn to the quantization of a massless scalar field in the presence of 
collapsing matter. The goal is to compute the flux of particles on J + given 
some initial state defined on J~ . To completely determine this flux on all of 
J~ + , one would have to solve the wave equation mode by mode in the region 
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bounded by J~ and E, where E is a constant t hypersurface which crosses 
the horizon to the future of the collapsing matter. This is clearly intractable, 
since the geometry inside the collapsing matter may be very complicated, 
and the scalar field might interact with the matter in an arbitrarily complex 
fashion. Fortunately, all of the dependence of the particle flux on J + on 
these factors dies out at sufficiently late times, and the radiation becomes, at 
least to leading order, completely independent of the details of the collapse 
process. Thus we shift our goal to calculating this late-time radiation and to 
showing that it is indeed universal. 

The first important observation is that the geometry is entirely smooth as 
long as one stays sufficiently far from the singularity at r = 0. In particular, 
for a black hole much larger than m p , the geometry is smooth near the horizon 
- in fact, the curvature can be made arbitrarily small by making the black 
hole arbitrarily massive. This strongly suggests the conlusion that however 
complicated the state of the field is after propagating through the matter, 
it should certainly appear nonsingular to inertial observers near the horizon. 
Of course, this presupposes that the initial state on J~ is nonsingular. 

Next, we examine the outgoing null geodesies in the region exterior to 
the collapsing matter. These obey: 



u = t- 2\j2Mr - r - 4Mlog (y / r/2M - 1) = constant. (3.21) 

Let us consider two geodesies, labelled by u\ and u 2 , which are separated 
in time by Am = u 2 — u\. On a constant t surface, E, the geodesies have a 
radial separation given by 

Au = n - r 2 + 2( v /2Mr 1 - j2Mr 2 ) + AM log ( ^ - V J^ \ , (3.22) 

\ y /r 2 ~-\ / 2mJ 

For large ui, u 2 , which corresponds to geodesies which reach J + at late 
times, the radial separation of the geodesies near the horizon is determined 
by 

Au^AMlogt ^-^ ) (3.23) 



or 

r 2 -2M ^ (n - 2M)e" Au/4M . (3.24) 

In other words, the outoging geodesies pile up along the horizon. Alterna- 
tively, we can note that all the geodesies which reach J + after u — u, where 
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u ^> M, were contained on the surface t = in the region between r = 2M 
and r = 2M + 4Me _M//4M . This has two important consequences. First, we 
see that the late time radiation is entirely determined by the state of the 
field at distances arbitrarily close to the horizon. Second, an outgoing wave 
suffers an arbitrarily large redshift in escaping from near the horizon to J + . 
So, putting it all together, we see that to compute the late-time radiation we 
need only consider nonsingular states in a region near the horizon. 

It is easy to construct a nonsingular state if we first define a new time 
coordinate. The trouble with the coordinate t used in ( p.3|) is that its flow is 
spacelike inside the horizon. However, this can be rectified if we define a new 



time, r, as the value of t along the curve dr + ^j2M/rdt = 0. r has the virtue 
of being nonsingular and timelike. We can choose a set of modes which are 
positive frequency with respect to r, and expand the field in terms of them. 
Then, by the equivalence principle and standard quantum field theory in 
flat space, the corresponding vacuum will have a nonsingular stress-energy 
as seen by inertial observers. Of course, this choice of state is not unique 
since other nonsingular states can be obtained by applying particle creation 
operators to this state. This is irrelevant if one is only interested in late-time 
radiation, since the excitations above the vacuum will eventually redshift 
away. 

Now at spatial infinity (more accurately: conformal infinity X + ) the vac- 
uum state is defined locally by the requirement that modes having positive 
frequency with respect to the variable u = t s — r* are unoccupied, where t s 
is Schwarzschild time and r* = r + 2Mln(r — 2M) is the tortoise coordinate. 
We wish to find the relationship between this requirement and the preceding 
one. 

The relationship between t and t s is 






(3.25) 



so that 



u = t e -r* = t- 2v / 2Mr - r - 4M ln(y/r - V2M) 



(3.26) 



and thus one finds that along a curve with dr + 




dt = 0, 




(3.27) 
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Because the last term on the right hand side is singular, the two definitions of 
positive frequency - with respect to u or to t - do not coincide. To remove the 
singularity, note that along any of the curves of interest e~ u / 4M has a simple 
zero at r = 2M, but is otherwise positive. Clearly then demanding positive 
frequency with respect to t along such curves requires positive frequency not 
with respect to u but rather with respect to 

U = -e""/ 4M . (3.28) 



In this way we have arrived at the famous Unruh boundary conditions 

To summarize, the appropriate construction near the horizon is to expand 
the field as 

<P = J [b^ u + ble- u ] (3.29) 

and take the state resulting from collapse to be \0u), where b u \0u) = 0. For 
simplicity, we have suppressed the field's dependence on the ingoing modes. 
It remains to describe this state in terms of particles defined by the modes 
e~ tuJu . Since u — > oo at the horizon, we need another set of modes in which 
to expand the field inside the horizon. These will be referred to as u™(u); 
their explicit form will not be needed. We then have 

' L: [(a w e'^ + ate^) Q(-U) + (a>» + Ckr) ©(^)]- (3-30) 



The relation between the modes for U < is given by 

'duo 



so 



yt/(M) = /|[«^ e "™ + ^ e H 5 U<0 ( 3 - 31 ) 

/OO 

due i{u,u{u)-uu) f (3_ 32 ) 
-oo 



or, since U(u) = — e 



-u/AM 



Aniuj .—Aniw 

otu'u = 4M e 2 * Mw T(l-4iMu) ; ^ u = 4M e -2nMu T {l+AiMu). 

uj' Uj' 

(3.33) 
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So 



2 _(4M) 3 7iuj 1 2 _{AMfnuj e" 87rMw 

|CVu;| " ^ 1 _ e -87rMw ' — X _ e -8wMw " 

(3.34) 

If we are only interested in the radiation which flows to infinity it is appro- 
priate to form a density matrix by tracing over the states inside the horizon. 
To find the average number of particles radiated we should integrate |/3 W / W | 2 . 
But because the black hole radiates for an infinite amount of time at a con- 
stant rate (in the current approximation), this yields infinity. To find the 
rate of emission, we can place the hole in a large box and use the density of 
states du/2-n for outgoing particles. For normalized modes the completeness 
relation, (|3.18|) , and the relation 

K>u,\ 2 = e 87rM 1/W (3.35) 

imply 

W|/Uf = e87rM l_ r (3.36) 
The rate of emission of particles in the range u to u + duo is then: 

F( u ) = — -— ^ • (3-37) 

This is precisely the rate of emission from a black body, in one dimension, at 
temperature T = l/8nM. F does not quite give the flux seen at infinity since 
some fraction, 1 — T(u), of the particles will be reflected by the spacetime 
curvature back into the hole. Thus, for the flux at infinity we write, 

FooM = -^P— (3.38) 

The preceding analysis shows that the average flux is that of a thermal 
body, not that the full density matrix is thermal. That the density matrix 
is exactly thermal is, in fact, easy to show by utilizing a clever trick due to 
Unruh [p6|] . We have not chosen to use this method here, and instead have 
gone through the rather laborious procedure of computing the Bogoliubov 
coefficients directly, because only the latter method can be used when one 
wants to find corrections to the spectrum. 
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3.3.1 Reissner-Nordstrom 

A similar analysis can be carried out for the Reissner-Nordstrom geometry, 
the metric for which is, in the L = 1, R = r gauge: 



ds 1 = -dt 2 + (dr + ^2M/r - Q 2 /r 2 dt) 2 + r 2 {d9 2 + sm 2 6d<j) 2 ). (3.39) 

The geometry has two horizons which are located at the two values of r for 
which t goes from being timelike to spacelike, and vice versa. The outer 
horizon radius is 



R+ = M + \j M 2 - Q 2 (3.40) 
and the inner horizon radius is 



R_ = M - yM 2 - Q 2 . (3.41) 

It should be noted that the coordinates in ( |3.39| ) fail to cover the region 
inside the inner horizon. This will not pose any obstacle to determining the 
radiance, though, as the calculation only depends in the geometry in the 
vicinity of the outer horizon. 

In the present case, the state resulting from collapse is specified by de- 
manding the absence of particles positive in frequency with respect to t along 
the curve dr + ^jlMjr — Q 2 /r 2 . Now the calculation can proceed in exact 
analogy to the uncharged case. The resulting temperature is: 

T(MQ) 1 „ (3 42) 

The temperature vanishes for the extremal black hole M = Q. If Q > M, 
consideration of the global geometry reveals that there is no black hole at all, 
but rather a naked singularity, which we will refer to as the meta-extremal 
case. Finally, if the radiated particles are themselves charged, the factor 
governing the emission probability is not the Boltzmann factor, e - "^, but 

/ u -Qq/R+(M,Q) \ 
CXP { T(M,Q) ) < 3 ' 43 > 

where the particles carry charge q. The factor Q/R+ is the electrostatic 
potential at the outer horizon, and appears as a chemical potential causing 
particles with the same sign of charge as the hole to be preferentially emitted. 
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3.3.2 Breakdown of the Semiclassical Approximation 

In the preceding treatment we have followed Hawking and considered the 
propagation of free fields on a classical background geometry, ignoring the 
back reaction altogether. It is usually assumed that the dominant effects of 
the back reaction can be incorporated by allowing the black hole to lose mass 
quasistatically, so that the radiation at infinity is thermal but with a slowly 
varying temperature. Indeed, this is how a normal thermal body is expected 
to behave during cooling. While this scenario seems quite plausible for a 
wide range of configurations, both for black holes and normal objects, there 
are notable instances where it is certain to be invalid, even qualitatively |27] . 
In particular, if the emission of a single typical quantum induces a relatively 
large change in temperature, a quasistatic description is clearly inappropri- 
ate. Since a typical quantum has energy equal to T, the temperature of the 
hole, we expect non-trivial back reaction effects to become important when 

T|£-T. (3.44) 

For a Schwarzschild black hole this happens when M ~ m p , whereas for 
a Reissner- Nordstrom hole it happens when M m Q. In the Schwarzschild 
case, the breakdown occurs in a regime of high curvatures and large quantum 
fluctuations, presumably inaccessible to semiclassical methods. Not so, how- 
ever, for a large near-extremal Reissner-Nordstrom hole, M ~ Q ^> m p , for 
which the curvature near the horizon remains small, suggesting that mean- 
ingful calculations can be performed without the need for a full theory of 
quantum gravity. In the next chapter we perform such a calculation and see 
explicitly that the radiance from the near extremal hole is markedly different 
from what free field calculations indicate. 



3.4 Effect of Fluctuations 

The semiclassical approach to back reaction effects in quantum gravity is to 
compute < T^ > for a field quantized on a classical background geometry, 
and then to use this quantity on the right hand side of Einstein's equations, 

- -g^Tl = 8tt < T^ > . (3.45) 
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The (generally intractable) problem then becomes to compute < T^ v > for 
an arbitrary metric. One expects this approach to provide a reasonable 
description of the gross effects of back reaction provided that the fluctuations 
in T^ v are not too large. For example, replacing < T^ v > by 



< T 2 

^- ± [US 



> 

should not lead to any substantial change in the geometry if ( |3.45|) is to be 
physically relevant. 

In order to be quantitative, we shall consider the computation of stress 
energy fluctuations in the moving mirror model, which is discussed in ^9] , 
|30| . We consider this model here because it provides a vivid demonstration 
of just how misleading the semiclassical approximation can be, and because 
we will referring to it later in another context. By "moving mirror model", 
we simply mean quantizing a scalar field in two dimensions subject to the 
condition that the field vanish along the timelike worldline z m (t). We will 
use the coordinates 

u = t- z ; v = t + z. (3.46) 
Without the mirror present the propagator is 



G(l, 2) = (O|T[0(1)0(2)]|O) = — log [{u 2 - Ul )(v 2 - Vl )} (3.47) 



and for a mirror at rest, z m (t) = 0: 



67(1,2) = — log 



(u 2 - toiXga - vi] 
(u 2 - vi){v 2 - u x ] 



(3.48) 



Now consider a general mirror trajectory, which we shall write in terms of 
u, v as v m (u). It is easiest to proceed by defining new coordinates in terms of 
which the mirror is at rest. Working in the region to the right of the mirror 
and defining 

U(u)=v m (u) ; V(v)=v, (3.49) 



we see that the trajectory v 



corresponds to U = V, which is the 



desired result. In the new coordinates the metric is 



(U'iu))- 1 dUdV, 



(3.50) 



and because of conformal invariance, the propagator has the same form as 
before: 



£7(1,2) 



1 

47T 



log 



(t/M-t/M)fo2-m; 

(U(u 2 ) - ui) (v 2 - v x ) 



(3.51) 



CHAPTER 3. QUANTUM FIELD THEORY NEAR BLACK HOLES 30 



The basic object we will be considering is the renormalized propagator de- 
fined by subtracting ( [3.471 ) from (|3.51|) : 



G*(l,2) = -loj 



U{u 2 ) - U{u x ) 



(u 2 - Ui)(U(u 2 ) - vi)(v 2 -U(ux)) 



(3.52) 



This prescription simply corresponds to normal ordering with respect to stan- 
dard Minkowski space creation and annihilation operators. < T^ u > can now 
be found by operating on Gr(1, 2) according to 



so 



< >-- 



lim 



d d 



d d 
np ^(1,2). 



In the present case the only nonzero element is 



d d 
lim — — — — Gr(1, 2) 



1^2 dui du 2 
which, after some algebra, is found to be 



(3.53) 
(3.54) 

(3.55) 



< T >- 

^ Mil 



l 2 J Vm ^ u) du^v m {uY 



(3.56) 



The result can be written in terms of the mirror trajectory z m (t) using 

dv m (u) 1 + z m 



yielding 



^ Tun 



du 



1 (1 



1 - z' 



d 



'\-z- 



(3.57) 



(3.58) 



12vr (1 - z m ) 2 dt 

These expressions hold in the region to the right of the mirror. The expres- 
sions to the left of the mirror are obtained by interchanging u and v. 

These results can be used to find the radiation reaction force on the 
mirror. The force four-vector is, by energy conservation: 



(3.59) 
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where AT^ = T fiu (+) — T^ v {—) is the difference in T^ evaluated on the right 
and left hand sides of the mirror, and v v is the mirror's velocity four-vector. 
We obtain: 

F^'^-^A) (3.60) 
12tt \ dr 2 dr 2 J y J 

where r is proper time. It is amusing to compare this to the radiation reaction 
force on a point charge in classical electrodynamics, 

This expression, containing as it does third derivatives of position of with 
respect to time, has some well known peculiar features such as runaway 
solutions and pre-acceleration. These carry over to the present case, although 
we hasten to add that a realistic mirror has a high frequency cutoff above 
which it becomes transparent; the worried reader need not be concerned with 
the hazards of runaway mirrors. 

< T^ u T a/ 3 > is found by operating on the renormalized four point function, 
which by Wick's theorem is: 

G R (1, 2, 3, 4) = G R (1, 2)G R (3A) + G R (1, 3)G R (2, 4) + G R (1, 4)G R (2, 3). 

(3.62) 

We find, in the region to the right of the mirror: 

< {T uu f >= 3 < T uu > 2 ; < (T vv ) 2 >= 

1 (dv m /du) 2 

< J-uu-Lvv >— q o~I 7 M4 (3.63) 

The fluctuations become infinite at the position of the mirror, v = v m (u). 
To see how this might affect the motion of the mirror we can consider the 
fluctuations in the radiation reaction force. Specifically, let us look at the 
force normal to the mirror, 

F N = -FX (3.64) 
where n^n^ = — 1 and n^v^ = 0. The mean squared value is found to be: 

< F 2 N >= n^n a v^[< T^(-)T a(3 (-) > + < T^(+)T^(+) > 
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- 2 < T aP {-) X T^(+) >] . (3.65) 
This is clearly infinite; for instance, the second term contributes 

f (dv m /duf 



n u n v v u v v < Tw (+)T W (+) > = 



87T 2 [^-^ m (^ 



(3.66) 



v=v m (u) 



The fluctuations diverge even for a stationary mirror. In the spirit of the 
semiclassical approach we might have attempted to model the back reaction 
by solving 

map = Fp. (3.67) 

However, in the case of the moving mirror this can not be interpreted as 
the leading term of a fully quantum mechanical treatment, since the fluctua- 
tions diverge. Instead, it seems more likely that a fully quantum mechanical 
treatment simply does not exist. 

While these effects are most important for the moving mirror, there is no 
reason to believe that fluctuations are so violent in the case of a black hole 
much larger then the Planck mass. The point is simply that spacetime is 
locally flat near the horizon on a scale much larger than the Planck length. 
The divergences in the mirror model arose from the sharp boundary condition 
at the mirror, a feature which has no analog in the black hole case. 



Chapter 4 

Self-Interaction Corrections 



Black hole radiance |3| was originally derived in an approximation where 
the background geometry was given, by calculating the response of quantum 
fields to this (collapse) geometry. As we have seen, in this approximation the 
radiation is thermal, and much has been made both of the supposed depth 
of this result and of the paradoxes that ensue if it is taken literally For if 
the radiation is accurately thermal there is no connection between what went 
into the hole and what comes out, a possibility which is difficult to reconcile 
with unitary evolution in quantum theory - or, more simply, with the idea 
that there are equations uniquely connecting the past with the future. To 
address such questions convincingly, one must go beyond the approximation 
of treating the geometry as given, and treat it too as a quantum variable. This 
is not easy, and as far as we know no concrete correction to the original result 
has previously been derived in spite of much effort over more than twenty 
years. Here we shall calculate what is plausibly the leading correction to the 
emission rate of single particles in the limit of large Schwarzschild holes, by 
a method that can be generalized in several directions, as we shall outline. 

There is a semi-trivial fact about the classic results for black hole radia- 
tion, that clearly prevents the radiation from being accurately thermal. This 
is the effect that the temperature of the hole depends upon its mass, so that 
in calculating the "thermal" emission rate one must know what mass of the 
black hole to use - but the mass is different, before and after the radiation! 
(Note that a rigorous identification of the temperature of a hot body from 
its radiation, can only be made for sufficiently high frequencies, such that 
the gray-body factors approach unity. But it is just in this limit that the 
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ambiguity mentioned above is most serious.) As we have emphasized, this 
problem is particularly quantitatively acute for near-extremal holes — it is a 
general problem for bodies with finite heat capacity, and in the near-extremal 
limit the heat capacity of the black hole vanishes. 

To resolve the above-mentioned ambiguity, one clearly must allow the 
geometry to fluctuate, namely to support black holes of different mass. An- 
other point of view is that one must take into account the self- gravitational 
interaction of the radiation. 

4.1 The Thin Shell Model 

To obtain a complete description of a self-gravitating particle it would be 
necessary to compute the action for an arbitrary motion of the particle and 
gravitational field. While writing down a formal expression for such an ob- 
ject is straightforward, it is of little use in solving a concrete problem due to 
the large number of degrees of freedom present. To arrive at a more work- 
able description of the particle-hole system, we will keep only those degrees 
of freedom which are most relevant to the problem of particle emission from 
regions of low curvature. The first important restriction is made by consider- 
ing only spherically symmetric field configurations, and treating the particle 
as a spherical shell. This is an interesting case since black hole radiation into 
a scalar field occurs primarily in the s-wave, and virtual transitions to higher 
partial wave configurations are formally suppressed by powers of 

Before launching into the detailed calculation, which becomes rather in- 
tricate, it seems appropriate briefly to describe its underlying logic. After 
the truncation to s-wave, the remaining dynamics describes a shell of matter 
interacting with a black hole of fixed mass and with itself. (The mass as 
seen from infinity is the total mass, including that from the shell variable, 
and is allowed to vary. One could equally well have chosen the total mass 
constant, and allowed the hole mass to vary.) There is effectively one degree 
of freedom, corresponding to the position of the shell, but to isolate it one 

1 Since we do not address the ultraviolet problems of quantum gravity these corrections 
are actually infinite, but one might anticipate that in gravity theory with satisfactory 

2 A 2 

ultraviolet behavior the virtual transitions will supply additive corrections of order u Mi , 
where A is the effective cutoff, and M is the mass of the black hole, but will not alter the 
exponential factors we compute. 
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must choose appropriate variables and solve constraints, since the original 
action superficially appears to contain much more than this. Having done 
that, one obtains an effective action for the true degree of freedom. This 
effective action is nonlocal, and its full quantization would require one to 
resolve factor-ordering ambiguities, which appears very difficult. Hence we 
quantize it semi-classically, essentially by using the WKB approximation. 
After doing this one arrives at a non-linear first order partial differential 
equation for the phase of the wave function. This differential equation may 
be solved by the method of characteristics. According to this method, one 
solves for the characteristics, specifies the function to be determined along 
a generic initial surface (intersecting the characteristics transversally) , and 
evolves the function away from the initial surface, by integrating the action 
along the characteristics. (For a nice brief account of this, see ||31|| .) 

When the background geometry is regarded as fixed the characteristics for 
particle motion are simply the geodesies in that geometry, and they are essen- 
tially independent of the particle's mass or energy — principle of equivalence 
- except that null geodesies are used for massless particles, and timelike 
geodesies for massive particles. Here we find that the characteristics depend 
on the mass and energy in a highly non-trivial way. Also the action along 
the characteristics, which would be zero for a massless particle and propor- 
tional to the length for a massive particle, is now a much more complicated 
expression. Nevertheless we can solve the equations, to obtain the proper 
modes for our problem. 

Having obtained the modes, the final step is to identify the state of the 
quantum field — that is, the occupation of the modes — appropriate to the 
physical conditions we wish to describe. We do this by demanding that a 
freely falling observer passing through the horizon see no singular behavior, 
and that positive frequency modes are unoccupied in the distant past. This, it 
has been argued, is plausibly the appropriate prescription for the state of the 
quantum field excited by collapse of matter into a black hole, at least in so far 
as it leads to late-time radiation. Using it, we obtain a mixture of positive- 
and negative- frequency modes at late times, which can be interpreted as 
a state of radiation from the hole. For massless scalar particles, we carry 
the explicit calculation far enough to identify the leading correction to the 
exponential dependence of the radiation intensity on frequency. 
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4.1.1 Effective Action 

We now derive the Hamiltonian effective action for a self-gravitating particle 
in the s-wave. First, we would like to explain why the Hamiltonian form of 
the action is particularly well suited to our problem. As explained above, 
our physical problem really contains just one degree of freedom, but the orig- 
inal action appears to contain several. The reason of course is that Einstein 
gravity is a theory with constraints and one should only include a subset of 
the spherically symmetric configurations in the physical description, namely 
those satisfying the constraints. In general, in eliminating constraints Hamil- 
tonian methods are more flexible than Lagrangian methods. This appears to 
be very much the case for our problem, as we now discuss. 

In terms of the variables appearing in the Lagrangian description, the 
constraints have the form 



= 0, 



where f is the shell radius, and ' represents 4. When applied to the spheri- 
cally symmetric, source free, solutions, one obtains the content of Birkhoff 's 
theorem - the unique solution is the Schwarzschild geometry with some mass, 
M. Since this must hold for the regions interior and exterior to the shell (with 
a different mass M for each), and since M must be time independent, we see 
that only those shell trajectories which are "energy conserving" are compat- 
ible with the constraints. This feature makes the transition to the quantum 
theory rather difficult, as one desires an expression for the action valid for an 
arbitrary shell trajectory. This defect is remedied in the Hamiltonian formu- 
lation, where the constraints are expressed in terms of momenta rather than 
time derivatives, 

Ch [r,P',9ij,^ij] = 0. 

At each time, the unique solution is again some slice of the Schwarzschild 
geometry, but the constraints no longer prevent M from being time depen- 
dent. Thus, an arbitrary shell trajectory f(t), pit), is perfectly consistent 
with the Hamiltonian form of the constraints, making quantization much 
more convenient. 

As before, we begin by writing the metric in ADM form: 



ds 2 = -N\t,rfdt 2 +L(t,rf[dr+N r (t,r)dt} 2 +R(t,rf[d9 2 +sm9 2 d(P 2 } (4.1) 
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In considering the above form, we have restricted ourselves to spherically 
symmetric geometries at the outset. With this choice of variables, the action 
for the shell is written 



S s = -m J ^-g^dxudx" = ~ m J dt\JN t2 - L 2 (r + iV r ) 2 , (4.2) 

m representing the rest mass of the shell, and the carets instructing one to 
evaluate quantities at the shell (g^ = ^(t, f)). 
The action for the gravity-shell system is then 

S = J d 4 x \J — g 1Z — m Jdt \j (N 1 ) 2 — L 2 (r + N r ) 2 + boundary terms 

* (4-3) 
and can be written in canonical form as 

S = Jdtpir+ Jdtdr[7i R R+7i L L-N t (H s t +nf)-N r {H s r +nf)}-JdtM ADM 

(4-4) 

with 



H s t = V (p/L) 2 + m 2 5{r - r) ; H s r = -p S(r - f) (4.5) 
G W ttlttr fRRf\ R' 2 L G 

where ' represents j-. Madm is the ADM mass of the system, and is numer- 
ically equal to the total mass of the combined gravity-shell system. 

We now wish to eliminate the gravitational degrees of freedom in order to 
obtain an effective action which depends only on the shell variables. To ac- 
complish this, we first identify the constraints which are obtained by varying 
with respect to N* and iV r : 

H t = H s t +H? = ; H r = H s r + H ( ; = 0. (4.7) 



By solving these constraints, and inserting the solutions back into ( [4.4|) we 
can eliminate the dependence on txr and 7Tl. We first consider the linear 
combination of constraints 

= %-Ut + ^rU r = -M' + ^Ht + ^rH S r (4.8) 
L RL L RL 
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where 

M = ^- + --^C- (4.9) 
2R 2 2L 2 v ' 

Away from the shell the solution of this constraint is simply A4 = constant. 

By considering a static slice {n L = tt r = 0), we see that the solution is a 

static slice of the Schwarzschild geometry with M. the corresponding mass 

parameter. The presence of the shell causes M. to be discontinuous at f , so 

we write 

M = M r < f 

M = M + r>r. (4.10) 

As there is no matter outside the shell we also have Madm — M+. Then, 
using ( [4.8| ) and ( |4.9| ) we can solve the constraints to find til and n R : 

tt l = R^(R'/L) 2 - 1 + 2M/R ; *r = ^/l r<t 

n L = Ry/(R'/L)* - 1 + 2M+/R ; ir R = ^ L ; r > f. (4.11) 

The relation between M + and M is found by solving the constraints at the 
position of the shell. This is done most easily by choosing coordinates such 
that L and R are continuous as one crosses the shell, and ttr^l are free of 
singularities there. Then, integration of the constraints across the shell yields 

ttl(^ + e) - n L {r - e) = -p/L 



R'{f + e) - R\f - e) = -i ^Jp 2 + m 2 L 2 (4.12) 

R 

Now, when the constraints are satisfied a variation of the action takes the 
form 

dS = pdf + J dr(n R 5R + n L 5L) - M + dt (4.13) 

where ttr^l are now understood to be given by ( |4.11| ), and M + is determined 
by solving (|4.12f ). We wish to integrate the expression ( 4.13 ) to find the 



action for an arbitrary shell trajectory. As discussed above, the geometry 
inside the shell is taken to be fixed (namely, M is held constant) while the 
geometry outside the shell will vary in order to satisfy the constraints. It is 
easiest to integrate the action by initially varying the geometry away from 
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the shell. We first consider starting from an arbitrary geometry and varying 
L until 7Tr = ttl = 0, while holding f,p, R, L fixed: 

JdS = J?° dr tf =0 8L n L 



J r r 2dr J^SL RJ{R'/LY - 1 + 2M/R + f™ e dr tf =0 5L RJ(R'/Lf - 1 + 2M+/R 



Ir~ e dr 

1 mm 



RLMR'/L) 2 - 1 + 2M/R + RR' log 



R' /L— y/(R' /L) 2 -1+2M/R 



^\1-2M/R\ 



RLJ(R'/L) 2 - 1 + 2M+/R + /^'log 



R'/L-y/(R'/L) 2 ~l+2M + /R 
y/\l-2M+/R\ 



(4.14) 

where the lower limit of integration, r m j n , properly extends to the collapsing 
matter forming the black hole; its precise value will not be important. We 
have discarded the constant arising from the lower limit of the L integration. 
In the next stage we can vary L and R, while keeping tt r l = 0, to some set 
geometry. Since the momenta vanish, there is no contribution to the action 
from this variation. 

It remains to consider nonzero variations at the shell. If an arbitrary 
variation of L and R is inserted into the final expression of ( 4.14 ) one finds 



dS 



dr [itrSR + tcl5L] — 



OS . dS / 
r + e) — (r — e) 



dR 



OR' 



dR+^^dM^ 



DM, 



(4.15) 



Since R' is discontinuous at the shell, 



OS ^ . OS 
(r + e) — (r 



dR' K ' dR n 
is nonvanishing and needs to be subtracted in order that the relations 



5S 



5S 



SR n ' SL 
will hold. From ( |4.14j ), the term to be subtracted is 
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'OS 
OR' 



(r + e) 



OS 
OR' 



r - e 



dR 



R'(f - e)/L - J(R'(r - e)/L) 2 - 1 + 2M/R 



+ dRR\og 



dRRlog 

R'(f + e)/L - J(R'(f + e)/L) 2 - 1 + 2M/R 



1 - 2M/R\ 



1 - 2M+/R\ 



(4.16) 



Similarly, arbitrary variations of L and R induce a variation of M + causing 
the appearance of the final term in ( 4.15 ). Thus we need to subtract 



dS 
dM7 



I r+e 1 - 2M+/R 4 



(4.17) 



Finally, we consider variations in p, f, and t. t variations simply give dS = 
—M + dt. We do not need to separately consider variations of p and f, since 
when the constraints are satisfied their variations are already accounted for 
in our expression for S, as will be shown. 

Collecting all of these terms, our final expression for the action reads 



S 



dr 



RLMR'/L) 2 - 1 + 2M/R + RR' log 



R'/L - J(R'/L) 2 - 1 + 2M/R 



1 - 2M/R\ 



+ / dr 

Jf+e 



RLJ{R'/L) 2 - 1 + 2M+/R + RR' log 



R'/L - J (R'/L) 2 - 1 + 2M+/R 



1 - 2M + /R\ 



dt RR 



+ log 



R'{f - e)/L - J(R'(f - e)/L) 2 - 1 + 2M/R 



1 - 2M/R\ 



log 



R'(f + e)/L - J(R'(r + e)/L) 2 - 1 + 2M/R 



1 - 2M + /R\ 
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+ dt 



- LyJ(R'/L) 2 -l + 2M + /R 
dr — ; M 4 

f+e 1 - 2M+/R 4 



dtM, 



(4.18) 



To show that this is the correct expression we can differentiate it; then it can 
be seen explicitly that when the constraints are satisfied ( |4.13| ) holds. 

We now wish to write the action in a more conventional form as the 
time integral of a Lagrangian. As it stands, the action in ( |4.18[ ) is given for 
an arbitrary choice of L and R consistent with the constraints. There is, 
of course, an enormous amount of redundant information contained in this 
description, since many L's and i?'s are equivalent to each other through a 
change of coordinates. To obtain an action which only depends on the truly 
physical variables p, f we make a specific choice for L and R, ie. choose a 
gauge. In so doing, we must respect the condition 



R'(f + e) - R'(f - e) 



— -VP 2 + m 2 L 2 
R 



which constrains the form of R' arbitrarily near the shell. Suppose we choose 
R for all r > f; then R'(f — e) is fixed by the constraint, but we can still 
choose R for r < f — e, in other words, away from the shell. We will let R' K 
denote the value of R' close to the shell but far enough away such that R is 
still freely specifiable. We employ the analogous definition for R' >: except in 
this case we are free to choose R> = R'(f + e). 

In terms of this notation the time derivative of S is 



L- — 
~ ~dt 



rRL 



(R'JL) 2 - 1 + 2M/R - x/iR^/L) 2 - 1 + 2M+/R 



-RRhg 



R'(f - e)/L - J{R'{f - e)/L) 2 - 1 + 2M/R 



+ 



R'JL 
dr[irjiR + ttlL] - 



{RJLf 



2M/R 



dr[n R R + n L L} - M + . 



(4.19) 



At this point we will, for simplicity, specialize to a massless particle (m = 0) 
and define r\ — ± = sgn(p). Then the constraints (|4.12|) read 
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(R'(f - e)/L) 2 - 1 + 2M/R = \/(R'(r + e)/L) 2 - 1 + 2M+/R + 



V 

LR 
(4.20) 



These relations can be inserted into ( j4.19|) to yield 
L = rRL 



(R^/L) 2 - 1 + 2M/R - y (R^/L) 2 - 1 + 2M+/R 



—r]RR log 



R'JL - W^/L) 2 - 1 + 2M + /R 



R'JL - r ] J(R' < /L) 2 - 1 + 2M/R 



+ / dr[7T R R + n L L]+ dr[ir R R + ir L L]- M 



(4.21) 



Now we can use the freedom to choose a gauge to make (|4.21|) appear as 
simple as possible. It is clearly advantageous to choose L and R to be time 
independent, so krR + tclL = 0. Also, having R' = L simplifies the ex- 
pressions further. Finally, it is crucial that the metric be free of coordinate 
singularities. A gauge which conveniently accommodates these features is 

L = 1 ; R = r 

The L = 1, R = r gauge reduces the Lagrangian to 



L = r[V2Mf — J2M + f] — r/rr log 



■t - T]V2M 

where M + is now found from the constraints ( 4.2(J| ) to be related to p by 

M . - M 



(4.22) 



r] - J2M + /r 



The canonical momentum conjugate to f obtained from 4.22 



is 



Pc 



dL i 

— r- = V2Mf - J2M + f - rjf lof 
df 



>r - r]V2M 

in terms of which we write the action in canonical form as 

S = 



dt[p c f - M 4 



(4.23) 



(4.24) 



(4.25) 



which identifies M + as the Hamiltonian. We should point out that M + is the 
Hamiltonian only for a restricted set of gauges. If we look back at (|4.21| ) we 
see that the terms 7TrR + tc^L will in general contribute to the Hamiltonian. 
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4.1.2 Quantization 



In this section we discuss the quantization of the effective action ( 4.25|) . First, 



it is convenient to rewrite the action in a form which explicitly separates out 
the contribution from the particle. We write 

M + = M- Pt 

so r 

S= dt\p c r+ Pt \ (4.26) 



and the same substitution is understood to be made in (|4.24|) . We have 



omitted a term, / dtM, which simply contributes an overall constant to our 
formulas. In order to place our results in perspective, it is useful to step 
back and consider the analogous expressions in flat space. Our results are an 
extension of 



p = ±^p t 2 - m 2 (4.27) 
S = Jdt[pr+p t ]. (4.28) 

Indeed, the G — > limit of Q4.24j) , (|4.25|) yields precisely these expressions 
(with 

m — 0). To quantize, one is tempted to insert the substitutions p — > —i-Jr, 



p t — > —i§t into ( 4.27p , so as t° satisfy the canonical commutation relations. 



This results in a rather unwieldy, nonlocal differential equation. In this trivial 
case we know, of course, that the correct description of the particle is obtained 
by demanding locality and squaring both sides of (|4.27|) before substituting 



p and pt- So for this example it is straightforward to move from the point 
particle description to the field theory description, i.e. the Klein-Gordon 
equation. Now, returning to ( |4.24| ) we are again met with the question of 
how to implement the substitutions p — > —id. In this case the difficulty is 
more severe; we no longer have locality as a guiding criterion instructing us 
how to manipulate Q4.24 ) before turning the p's into differential operators. 



This is because we expect the effective action ( |4.26| ) to be nonlocal on phys- 
ical grounds, as it was obtained by including the gravitational field of the 
shell. 

There is, however, a class of solutions to the field equations for which this 
ambiguity is irrelevant to leading order, and which is sufficient to determine 
the late-time radiation from a black hole. These are the short-wavelength 
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solutions, which are accurately described by the geometrical optics, or WKB, 
approximation. Writing these solutions as 

0(t,r) = e iS{t ' r \ 

the condition determining the validity of the WKB approximation is that 

\dS\ » \d 2 S\V 2 , \&s\ 1 / 3 ... 

and that the geometry is slowly varying compared to S. In this regime, 
derivatives acting on (p(t, r) simply bring down powers of dS, so we can 
make the replacements 

dS dS 

Pc -7T- j Pt TTT 

or Ot 

and obtain a Hamilton- Jacobi equation for S. Furthermore, it is well known 
that the solution of the Hamilton- Jacobi equation is just the classical action. 
So, if f(t) is a solution of the equations of motion found by extremizing ( |4.26| ), 
then 

S (t, f(t)) = S (0, f(0)) + fdt \p c (f(t)) 'fit) + pt] (4.29) 

J o L J 

where 

Pc (0,f) = — (0,f). (4.30) 



Since the Lagrangian in ( |4.26|) has no explicit time dependence, the Hamilto- 



nian p t is conserved. Using this fact, it is easy to verify that the trajectories, 
r(i), which extremize (|4.26| ) are simply the null geodesies of the metric 



ds z = -dt 2 +[dr + J^±dt I . (4.31) 



From ( p.8[ ) the geodesies are: 

ingoing: t + f(t) + 2J2M + f{t) + 4M+ log [Jr~if) + Jwi 



f(0) + 2 v /2M+r(0) + 4M + log [Jf(0) + J2M. 
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outgoing: t - f(t) - 2J2M + f(t) - 4M+ log [Jf(t) 



2M 4 



f (0) - 2 v /2M+f(0) - 4M+log [Jf(0) - J2M+]. (4.32) 



M + , in turn, is determined by the initial condition S(0, r) according to ( |4.24j ) 
and ( p0|) : 



ingoing: — (0,f(0)) = ^2Mr(0) - ^2M + f(0) + f(0) log 



: (0) + V2M 



outgoing: —(0, f(0)) = ^2Mf(0)-^2M+f(0)-f(0) log 



Finally, we can use this value of M + to determine p c (t): 



f (0) - 



yr(0) - V2M 
(4.33) 



ingoing: p c (t) = J2Mf(t) - J2M + f(t) + f (t) log 



f{t) + vOT; 



yr(t) + V2M 



f(t) - ^2Ml 



outgoing: pJt) = J2Mf(t) - J2M + f(t) - f(t) log 

!r(t) - V2M 

(4.34) 

These formulas are sufficient to compute S(t,r) given 5(0, r). 

As will be discussed in the next section, the relevant solutions needed to 
describe the state of the field following black hole formation are those with 
the initial condition 

S(0,r) = kr k>0 (4.35) 

near the horizon. Here, k must be large 1/M) if the solution is to be 
accurately described by the WKB approximation. In fact, the relevant fc's 
needed to calculate the radiation from the hole at late times become arbi- 
trarily large, due to the ever increasing redshift experienced by the emitted 
quanta as they escape to infinity. We also show in the next section that 
to compute the emission probability of a quantum of frequency u, we are 
required to find the solution for all times in the region between r = 2M 
and r = 2{M + iS). That said, we turn to the calculation of S(t,r) in this 



region, and with the initial condition ([4.35 ). The solutions are determined 
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from ( f4.29| ), ( |4.32| )-( p~34| ). Because of the large redshift, we only need to keep 
those terms in these relations which become singular near the horizon. We 
then have for the outgoing solutions: 



S(t,r) = kf(0) — / <if flog 

Jr(0) 



r - ^/2Ml 



(M+ - M)t (4.36) 



t - 4M+ log [Vr - ^2M+] = -4M+ log [y/r(0) 
k = — f (0) log 



2M: 



f (0) - v/27W; 



f (0) - \PM 



(4.37) 
(4.38) 



To complete the calculation, we need to invert fl4.37]) and ((4.38j) to find 
M + and f (0) in terms of t and r, and then insert these expressions into ( 4.36|) . 
One finds that to next to leading order, 

Uk/2M' _ iy-t/4M' 



2M + = V2M + (Vr - V2M) 



1 _)_ ( e fe/2M' _ l\ e -t/4M' 
e (k/2M'-t/AM') 



where 



M' = M+ V2M(v^ - V2M)- 



,(k/2M-t/4M) 



(4.39) 



(4.40) 



I _|_ e (k/2M-t/AM) ' 

Plugging these relations into (|4.36| ) and keeping only those terms which con- 
tribute to the late-time radiation, one finds after some tedious algebra, 



S{t,r) = -(2M 2 -r 2 /2)log 



3 (fc/2M'-t/4M') 



(4.41) 



4.1.3 Results 

We will now discuss the application of these results to the problem of black 
hole radiance. The procedure is a slight modification of the one we discussed 
in the free field theory case. As before, the point is that there are two 
inequivalent sets of modes which need to be considered: those which are 
natural from the standpoint of an observer making measurements far from 
the black hole, and those which are natural from the standpoint of an observer 
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freely falling through the horizon subsequent to the collapse of the infalling 
matter. The appropriate modes for the observer at infinity are those which 
are positive frequency with respect to the Killing time, t. Writing these 
modes as 



<f>(t,r) reads 



4>(t,r) = J 



dk 



u k {r)e 



a^(r)e- ! ^ + 4<(r)e^ 



(4.42) 



These modes are singular at the horizon, 

du k 



dr 



oo as r 



2M. 



Symptomatic of this is that the freely falling observer would measure an 
infinite energy-momentum density in the corresponding vacuum state, 



(0 t |T^|0 t ) 



oo as 



2M 



where dk\O t ) = 0. However, we do not expect this to be the state resulting 
from collapse, since the freely falling observer is not expected to encounter 
any pathologies in crossing the horizon, where the local geometry is entirely 
nonsingular for a large black hole. To describe the state resulting from col- 
lapse, it is more appropriate to use modes which extend smoothly through 
the horizon, and which are positive frequency with respect to the freely falling 
observer. Denoting a complete set of such modes by Vk(t,r), we write 



4>(t,r) Jdk [b k v k (t,r) + b{v* k (t,r) 



(4.43) 



Then, the state determined by 



|0„) = 



results in a non-singular energy-momentum density at the horizon, and so is 
a viable candidate. The operators d k and b k are related by the Bogoliubov 
coefficients, 

d k = dk' \a kk >b k > + f3 kk >bl,] (4.44) 
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where 

1 



Pw = I dt e^vl, (t, r) . (4.45) 



2nu k {r) 

Note that here we compute the Bogoliubov coefficients by performing a t 
integration, rather than an integration over a spatial coordinate, as is con- 
ventional. We are forced to do this in the present case since we do not know 
the spatial dependence of the definite energy modes. The flux seen at infinity 
is 

Foo (u k ) = — K . 4.46 

2tt \akk'/Pkk'\ ~ 1 

Next, we consider the issue of determining the modes Vk(t,r). As stated 
above, we require these modes to be nonsingular at the horizon. Since the 
metric near the horizon is a smooth function of t and r, a set of such modes 
can be defined by taking their behaviour on a constant time surface, say 
t — 0, to be 

v k (0,r)^e ikr as r -»• 2M. 

This is, of course, the initial condition given in ( J4.35[ ) . Now, the integrals 
in ([4.45 ) determining the Bogoliubov coefficients depend on the values of 



Vk(t,r) at constant r. Since v k is evaluated in the WKB approximation, 
the highest accuracy will be obtained when r is as close to the horizon as 
possible, since that is where v k s wavelength is short. On the other hand, in 
calculating the emission of a particle of energy u>k, we cannot take r to be 
less than 2(M + Uk), since the solution Uk{r)e~ luJkt cannot be extended past 
that point. Therefore, we calculate the integrals with r = 2(M + Uk)- 

The results of the previous section give us an explicit expression for t> fc . 
From ( pETD , 



V k (t, 2(M + UJ k )) = e WW+*k)) = e .(4M^ + 2^)lo g [l + e(^/^'-/^')] (44?) 

where M' is 



{k/2M-t/AM) e {k/2M-t/AM) 



M> = M+V2M(J2(M + co k )-y/2M) , , , k/9M _ fUM , « M+u k 



X _|_ e (k/2M-t/4M) ' K l _|_ e {k/2M-t/AM) ' 

" (4.48) 



CHAPTER 4. SELF-INTERACTION CORRECTIONS 



49 



Then, the integrals are, 

dte t Wfc t e ±i(4Af Wfc ,+2 W » / )log[l+e(*'/ aM '-*/* M ')] j ( 4 _ 49 ) 



the upper sign corresponding to akk', and the lower to fikk'- We can compute 
the integrals using the saddle point approximation. It is readily seen that for 
the upper sign, the saddle point is reached when 

e (fc'/2M'-t/4M') 

so t is on the real axis. For the lower sign, the saddle point is 

e (fc'/2M'-t/4M') ^ -1/2, 

which, to zeroth order in u^, gives 

t = AinM + real 

and to first order in Uk, gives 

t = 4in(M - u k ) + real. 

Inserting these values of the saddle point into the integrands gives for the 
Bogoliubov coefficients, 



Olkk' 



e ^(M-w k )u) k ^ (4.50) 



Pkk' 

The flux of radiation from the black hole is given by (|4.46p , 

FooM = ^ = J M . (4.51) 

There is an alternative way of viewing the saddle point calculation, which 

provides additional insight into the physical origin of the radiation. Let us 
rewrite the integral ( |4.49|) as 

3 dt e iu >kt±iS{t,2(M+u k ,)) ^ (4.52) 
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The saddle point is given by that value of t for which the derivative of the 
expression in the exponent vanishes: 

dS 

u k ±—{t,2{M + u k )) = 0. 
at 

But dS/dt is just the negative of the Hamiltonian, 

dS 



dt 



p t = M-M+ 



so the saddle point equation becomes 

M + = M± u k . 



To find the corresponding values of t, we insert this relation into ([4.37 ) and 

(EHD : 



t = 4{M± u k ) log 



2(M + uo k ) + e - J2(M + oo k 



(4.53) 



-f(0) k>£ 



r(0)- J2{M±uj h ) 



: (0) - V2M 



(4.54) 



where we have written f = 2(M + u k ) + e to make explicit that f must lie 
outside the point where the solutions u k (r) break down. We desire to solve 
for t as k — > oo. For the upper choice of sign, we find from ( 4.54j) that 



V^(0) = y/2(M + u; k ) + 0(e~ fc / 2M ) 



which, from (|4.53|) , then shows that the corresponding value of t is purely 
real. 

For the lower choice of sign we have, 



y/m = p{M-u k ) - O (e" fc / 2M ) . 
Continuing t into the upper half plane, we find from ( 4.53 ) that 

t = Am (M - u k ) + real. 
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These results of course agree with our previous findings. 

The preceding derivation invites us to interpret the radiation as being 
due to negative energy particles propagating in imaginary time. The particles 
originate from just inside the horizon, and cross to the outside in an imaginary 
time interval 47r(M — Uk). This, perhaps, helps clarify the analogy between 
black hole radiance and pair production in an electric field, which, in an 



instanton approach p2| , is also calculated by considering particle trajectories 
in imaginary time. 

Finally, let us return to the question of thermality. One might have 
guessed that the correct exponential suppression factor could be the Boltz- 
mann factor for nominal temperature corresponding to the mass of the hole 
before the radiation, after the radiation, or somewhere in between. Thus one 
might have guessed that the exponential suppression of the radiance could 
take the form e _a;//Tbcforc , e _a; / Taftor , or something in between. Our result, 
to lowest order, corresponds to the nominal temperature for emission being 
equal to T after . 



4.2 Corrections to Charged Black Hole Ra- 
diance 

In this section, two additional things are done. First, we extend the calcu- 
lations to include a charged black hole, and charged matter. Although this 
step does not present any significant formal difficulties, the physical results 
we obtain are considerably richer than what we found in our previous calcu- 
lations involving neutral holes and shells. In the neutral case the final result 
could be summarized as a simple replacement of the nominal temperature 
governing the radiation by the Hawking temperature for the mass after radi- 
ation, so that the "Boltzmann factor" governing emission of energy u from 
a hole of mass M became 

g-w/Teff. _ e -u8Tr(M-u) ^ (4.55) 

Note that the argument of the exponential is not simply proportional to the 
energy u, so that the spectrum is not, strictly speaking, thermal. While 
the deviation from thermality is important in principle its structure, in this 
case, is rather trivial, and one is left wondering whether that is a general 
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result. Fortunately we find that for charged holes the final results are much 
more complex. We say "fortunately" , not only because this relieves us of the 
nagging fear that we have done a simple calculation in a complicated way, 
but also for more physical reasons. For one knows on general grounds that 
the thermal description of black hole radiance breaks down completely for 
near-extremal holes. One might anticipate, therefore, that something more 
drastic than a simple modification of the nominal temperature will occur - as 
indeed we find. A particularly gratifying consequence of the accurate formula 
is a form of "quantum cosmic censorship". Whereas a literal application of 
the conventional thermal formulas for radiation yields a non-zero amplitude 
for radiation past extremality - that is, radiation leaving behind a hole with 
larger charge than mass - we find (within our approximations) vanishing 
amplitude for such processes. 

Second, we discuss in a more detailed fashion the relationship between 
our method of calculation, which proceeds by reduction to an effective par- 
ticle theory, and more familiar approximations. We show that it amounts 
to saturation of the functional integral of the underlying s-wave field theory 
with one-particle intermediate states, or alternatively to neglect of vacuum 
polarization. It is therefore closely related to conventional eikonal approxima- 
tions. We demonstrate the reduction of the field theory to a particle theory 
explicitly in the related problem of particle creation by a strong spherically 
symmetric charge source, which is a problem of independent interest. 



4.2.1 Self- Interact ion Correction 

Our system consists of a matter shell of rest mass m and charge q interacting 
with the electromagnetic and gravitational fields. The corresponding action 
is 

S = f [-mJ-g. u dx^dx v + qAJx»\ + — f d^x^g [TZ - F^F^] (4.56) 

The gravitational contribution to the Hamiltonian is the same as in 
and the shell and electromagnetic contributions are 



H\= V (p/L) 2 + m 2 - qA t 5(r - f) ; H s r = -p S(r - f) (4.57) 



/EM iV*L7!X 2 



nr = -^^-A^ T (4.58) 
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To arrive at this form we have chosen a gauge such that A t is the only non- 
vanishing component of A^. Of course, we set A r = only after computing 
the canonical momentum 7T^ r . 

Constraints are found by varying the action with respect to N t , N r , and 

At, 

nj i/.s I n_/G i njEM r\ . n_j i/s i n/G r\ 

rit = ri t + /r t + /r t = U , /T r = ri r + rty, = U 

vr^ + g5(r-f) = 0. (4.59) 
7T/j can be eliminated by forming the linear combination of constraints 

= f 7U + ^H r = —M' + *{7ii + Hf M ) + (4-60) 

where 

7Tr, R RR 

We see from the Gauss' law constraint that — txa t if) is the charge contained 
within a sphere of size r, so we define: Q(r) = — 7u r (r) 

Now, if the shell was absent A4(r) and Q(r) would be given by 

r°° R'(r)H EM (r) 
M(r) = M- dr K > T ,\ V) ; Q(r) = Q (4.62) 
Jr L(r) 

with M and Q being the mass and charge of the black hole as seen from 
infinity. In the gauge L — 1, R — r these become 

M(r) = M-Q 2 /2r ; Q(r) = Q. (4.63) 

With the shell present we retain the expression ( |4.63| ) for the region inside 
the shell, r < f, whereas outside the shell we write (with L = 1, R = r), 

M(r) = M + -(Q + q) 2 /2r ; Q(r) = Q + q (4.64) 

where M + and Q + q are the mass and charge of the hole-shell system as 
measured at infinity. 

By using the constraints we can determine tir, 7Tl, and an expression for 
M_|_, in terms of the shell variables. These relations can then be inserted in 
the action to give an effective action depending only on the shell variables. 
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The calculation for the present case runs precisely parallel to the uncharged 
case, resulting in 

S = Jdt 'r (^2Mr - Q 2 - pM+r - (Q + q) 2 ) 



-Tjff log 



f - -(Q + q) 2 /2f 



Vf- - ri^/M - Q 2 /2f 



(4.65) 



where 77 = sgn (p), and we have now specialized to a massless shell (m = 0). 
The canonical momentum is then 



p c = ^2Mf - Q 2 -pM + f - (Q + g) 2 -r?r log 



f - 77^/M+ - (Q + g) 2 /2f 



v 7 ? - f^M - Q 2 /2r 



(4.66) 

We need wish to find the short wavelength solutions which are accurately 
described by the WKB approximation. Writing these solutions as v(t,r) = 
e *S{t,r) w ith S rapidly varying, we can make the replacements 



Pc 



S(t, r) satisfies the Hamilton- Jacobi equation, and so is found by comput- 
ing classical action along classical trajectories. We first choose the initial 
conditions for S(t,r) at t — 0: 



S q k (0,r) = kr. 



(4.67) 



We have a appended a subscript and a superscript to denote the initial con- 
dition and charge of the solution. The corresponding classical trajectory has 
the initial condition p c = k at t = 0. Sl(t, r) is then given by 



SZ(t,r) = fcf(0) + / dfp c (f) - (M+ - M)t. 

'r(O) 



(4.68) 



To determine the radiance from the hole we will will only need to consider 
the behaviour of the solutions near the horizon. Furthermore, only the most 
rapidly varying part of the solutions will contribute to the late-time radiation. 
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With this in mind, we can write the momentum as (choosing 77 = 1 for an 
outgoing solution) 



p c (f) fa —f log 



r-R + (M + ,Q + q) 



(r-R+(M,Q))(r-R-(M,Q)) 
so that the initial condition becomes 

(f (0) - R + {M + , Q + g))(f (0) - i?_(M + , Q + q)) 



k 



-f (0) k>£ 



(f (0) — R+(M, Q))(f (0) — R-(M, Q)) 



(4.69) 



(4.70) 



Similarly, the classical trajectory emanating from f (0) is given by approxi- 
mately, 



R + (M + , Q + q) - R_(M+, Q + q) 
-R4M + ,Q + q) 2 log 



R+(M + ,Q + q) 2 log 
r-R4M+,Q + q) 



r-R + (M + ,Q + q) 



r(0)-R-(M + ,Q + q) 
These trajectories are in fact null geodesies of the metric 



ds 2 



-dt 2 + (dr + J2M + /r - Q 2 dt) A 



r(0)-R+(M + ,Q + q) 
(4.71) 

(4.72) 



The relations Q4.70Q and (|4.71|) allow us to determine M + and f(0) in terms 
of the other variables, so that after integrating (|4.68 ) we can obtain an ex- 
pression for Sl(t, r) as a function of k, t, and r. 

We can now write down an expression for the field operator: 



4(t,r)= [dk[d k vl(t,r) + b{v^(t,rn 



(4.73) 



The modes vl(t,r) are nonsingular at the horizon, and so the state of the 
field is taken to be the vacuum with respect to these modes: 



0,) = h \0 V ) = 0. 



Alternatively, we can consider modes which are positive frequency with re- 
spect to the Killing time t. We write these modes as u q k {r)e~ %ulkt where the 
u q k (r) are singular at the horizon, r = R + (M + Uk, Q + q). Then 



far) = [dk [c ft uj(r)e-*** + fc 9 (r)V^]. 



(4.74) 
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The two sets of operators are related by Bogoliubov coefficients, 

c k = jdk [a kk id k > + Pkk'b{>}- (4.75) 



From ( |4.73| , |4.74| ) a kk i and [3kk> are found to be 



f r°° 
2-Ku\(r) J-oo 

1 r°° 

Afc ' = 2^W dte^%,%r)*. (4.76) 



-co 



Here, r is taken to be slightly outside the horizon, r = R + (M + u k , Q + q) + 
e. These coefficients can be evaluated in the saddle point approximation. 
Recalling that vl(t,r) = e lS k( t > r \ the saddle point equation for a kk / becomes 

w* = = M%-M. (4.77) 

This leads to a purely real value of t for the saddle point. For f3 kk > we have 

u k = ^ r = M-M- + \ (4.78) 



From ( |4.70| , 4.71]) we find that the saddle point value for t has an imaginary 



part given by 

r /, s = 2/? + (M-u; fc ,Q-g) 2 ^ = 1 

m[s) R + (M-u k ,Q-q)-R_(M-u k ,Q-q) n 2T(M-u k ,Q-q) 

(4.79) 

Therefore, 

\Pkk'/a kk ,\ = i^^TT exp (u k /T(M -u k ,Q-q) + lm[Sj?(t a )*)) . (4.80) 

The terms in ST/, 9 which contribute to the second term in the exponent are 

/ drpJr) +Uklm(t 8 ). 
Jf(o) 
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Moo + sJM 2 -Q 2 (J(M - uj) 2 - (Q - g) 2 - ^M 2 - Q 2 



Using ( [4.69| - [171"1 ) this can be evaluated to give 

Ms;> q (t s y) 

resulting in 



2T(M-u,Q-q) R+(M,Q) 



(4.81) 



Pkk' 


I 




= exp 


Oikk 1 





VM' 2 - Q 2 [uj - sJ{M - uj) 2 - (Q - q) 2 + yJhP - Q 2 
T(M-u,Q-q) R+(M,Q) 

(4.82) 

This is the effective Boltzmann factor governing emission. Sufficiently far 
from extremality, when u, q \/M 2 — Q 2 , we can expand Q482D to eive 





1 


m exp 


O^kk' 





U) 



Qq 



+ 



M z q 2 +Q I uj I -2MQujq \ 



R+{M,Q) ^ 2(M 2 -Q 2 )R+(M,Q) \ 



T(M-u,Q-q) 
as compared to the free field theory result |§, 



(4.83) 



Pkk' 


I 




= exp 


Oikk' 





UJ — 



Qq 



R+(M,Q) 



(4.84) 



T(M,Q) 

Near extremality, the self-interaction corrections cause the emission to differ 



substantially from ( 4.84 ). 

We might ask whether it is possible to reach extremality after a finite 
number of emissions. Since T(M — u,Q — q) appears in the denominator of 
the exponent of (|4.82|) , the transition probability to the extremal state is in 
fact zero. We can also ask whether there are transitions to a meta-extremal 
(Q > M) hole. This would have rather dramatic implications as the meta- 
extremal hole is a naked singularity. To address this question we return 
to the saddle point equation ( |4.78|) . When Q > M, R + and R- become 
complex. From ( [4.70[ ) we see that a saddle point solution would require that 
k be complex, but we do not allow this since a complete family of initial 
conditions <Sfc(0, r) = kr was defined with k real. Therefore, in the saddle 
point approximation the extremal hole is stable. 

Modes with \f3/a\ > 1 formally require larger amplitudes for higher oc- 
cupation numbers, and thus require special interpretation. Considering for 
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simplicity the free field form of these coefficients, (|4.84 ), we see that such 
modes occur when uj < qQ/R +1 that is when the incremental energy gain 
from discharging the Coulomb field overbalances the cost of creating the 
charged particle. Under these conditions one has dielectric breakdown of the 
vacuum, just as for a uniform electric field in empty space. Since this physics 
is not our primary concern in the present work, we shall restrict ourselves 
to a few remarks. The occupation factor appearing in the formula for radia- 
tion in these "superradiant" modes is negative, but the reflection probability 
exceeds unity, so the radiation flux is positive as it should be. And in gen- 
eral the formulas for physical quantities will appear sensible, although Fock 
space occupation numbers are not. We can avoid superradiance altogether 
by considering a model with only massive charged fundamental particles, 
and holes with a charge/mass ratio small compared to the minimal value for 
fundamental quanta. 

Another interesting variant is to consider a magnetically charged hole 
interacting with neutral matter. In that case, one simply puts q = in the 
fomulae above (but Q ^ 0). One could also consider the interaction of dyonic 
holes with charged matter, and other variants (e.g. dilaton black holes) but 
we shall not do that here. 



4.3 Discussion 

We have arrived at our results by what may have appeared to be a some- 
what circuitous route. Inspired by a field theory question, we calculated the 
solutions of a single self-gravitating particle at the horizon, and then passed 
back to field theory by interpreting the solutions as the modes of a second 
quantized field operator. In this section we hope to clarify the logic of this 
procedure, and show that it is both correct and efficient, by demonstrating 
how a single particle action emerges from the truncation of a complete field 
theory. 

We can illustrate this explicitly if we consider the simpler model of spher- 
ically symmetric electromagnetic and charged scalar fields interacting in flat 
space. Our goal is to show that the propagator for the scalar field can be ex- 
pressed as a Hamiltonian path integral for a single charged shell. To achieve 
this, two important approximations will be made. The first is that the ef- 
fects of vacuum polarization will be assumed to be small, so we can ignore 
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scalar loop diagrams. The second is to assume that the dominant interactions 
involve soft photons, so that the difference in the scalar particle's energy be- 
fore and after emission or absorption of a photon is small compared to the 
energy itself. Thus we expect that our expression will be valid for cases 
where the scalar particle has a large energy, so that the energy transfer per 
photon is relatively small, and is far from the origin, so that the classical 
electromagnetic self energy of the particle is a slowly varying function of the 
radial coordinate. Field theory in this domain is in fact well described by the 
eikonal approximation, which implements the same approximations we have 
just outlined. What follows is then essentially a Hamiltonian version of the 
eikonal method. 

We start from the action 






(4.85) 



Defining the charge density 




(4.86) 



the solution of the Gauss' law constraint is 




(4.87) 



and so the scalar field Hamiltonian becomes 




(4.88) 



The fields are now written as second quantized operators: 





(4.89) 
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where u k = \/k 2 + m 2 , and we also have 0* = 0^ , tt^* = ix^ . To ensure that 
the field is nonsingular at the origin we impose the conditions d^ k = —d k , 

= — b k , and take the limits of all k integrals to be from — oo to oo. 

We now write the Hamiltonian in terms of the creation and annihilation 
operators. In doing so we shall normal order the operators, which corresponds 
to omitting vacuum polarization since we do not allow particle-antiparticle 
pairs to be created out of the vacuum. Also when evaluating 0' we shall use 
the geometrical optics approximation, (e lkr jr)' xs ike lkr /r, valid for k 3> 1/r. 
Then the quadratic part of the Hamiltonian becomes, 

f°°dr + r 2 0*'0' + m 2 0*0 =\ ( dku k [a\a k + b\b k \- (4.90) 

jo v 2 J 

Next we consider the interaction term. When evaluating this there will arise 



factors of yuj k //uj k . The essence of the soft photon approximation is that we 
replace these factors by 1, since we are assuming that Au/u 1 for the 
emission or absorption of a single photon. Then, after normal ordering, we 
can evaluate the charge density to be: 

p(r) = g/^[4^-^]e i(fc "^. (4.91) 

We now wish to calculate matrix elements of the Hamiltonian between 
one particle states. A basis of one particle states labelled by position is given 
by ' 

<lh e-^4|0>- (4.92) 



27T 

The free part of the Hamiltonian then has matrix elements 

(r 2 | l - J dku k [a\a k + 1%]}^) = u k [e ik ^~^ - e ik ^\ (4.93) 

The second term in the brackets corresponds to the path from r\ to r 2 which 
passes through the origin. These paths will not contribute to local processes 
far from the origin, so we drop this term. The matrix elements of the inter- 
action term for closely spaced points r\ and r 2 are: 
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Putting these together, we find the matrix elements of the Hamiltonian, 



Now we can follow the standard route which leads from matrix elements of 
the Hamiltonian to a path integral expression for the time evolution operator, 
with the result 



The action in the exponent is precisely that of a charged shell, with q 2 /2r 
being the electromagnetic self energy. 

We now discuss how this analysis can be applied to the case where we 
include gravitational interactions. The resulting field Hamiltonian is much 
more complex, and so we will not be able to explicitly calculate the effective 
shell action. However, the preceding derivation allows us to argue that were 
we to do so, we would simply derive the effective action obtained in section 
2. The nature of the black hole radiance calculation makes us believe that 
the approximations used to arrive at a shell action are justified. This is so 
because for a large (M ^> m p ) hole the relevant field configurations are short 
wavelength solutions moving in a region of relatively low curvature, and these 
are the conditions which we argued make the eikonal approximation valid. 

For simplicity, we will consider an uncharged self-gravitating scalar field. 
If we truncate to the s-wave we arrive at what is known as the BCMN model, 
originally considered in |T7J and corrected in [^|. The action is 




(4.95) 




(4.96) 




/ 



dt dr + ttrR + ti l L - JV*(Wf + H?) - N r (Hf + H G r ) 



I 



dtM J 



ADM 



(4.97) 



with 




(4.98) 



The analog of (|4.60| ) is now 




(4.99) 
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The Hamiltonian is 



H = Ma dm = M(oo) = M + / dr 



■Xt 







~L + R~L 



(4.100) 



To obtain an expression for H which depends only on and 7r^ we must 
choose a gauge and solve the constraints. We can obtain an explicit result if 
we choose the gauge R = r, n L = 0. Then, defining 

(4.101) 




L is determined from ( |4.99| ) 



so 

1 2M 



-2 



L 2 r 
which then leads to 



f Q r dr> h(r')/r> + 1 g -2 />' M^OA' T^/ g 2 J^dr" fc(r")/r" ( 4 _ 103) 

r Jo 



H = Me- 2 ^ drh{r)/r + rdrh{r)e- 2 l°° dr ' h{r ' )/r '. (4.104) 

Jo 



This generalizes the result of |26j to include a nonzero mass M for the pure 
gravity solution. To make a direct comparison with our work in the previous 
section, it would be preferable to obtain the Hamiltonian in L = 1,R = r 
gauge. This is more difficult and we do not know the explicit expression. For 
the moment, though, we are mainly interested in the qualitative structure of 



the Hamiltonian, and (|4.104|) will be sufficient for our purposes. The various 



nonlocal terms contained in the Hamiltonian ( |4.104j ) correspond to gravitons 
attaching onto the particle's worldline. If we expand the exponentials in 
( |4.104p , we see that there arise an infinite series of bi-local, tri-local, 



terms resulting from the non-linearity of gravity. Now we could, in principle, 
repeat the analysis which led to an effective shell action for the charged field 
in flat space. In that case the calculation could be done with only modest 
effort because there was only a single quartic interaction term. In the present 
case we would have to sum the infinite series of terms that arise; our point 
is that handling all of these terms is cumbersome, to say the least, and that 
it is much simpler to proceed as in section 2. 
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4.4 Multi-particle Correlations 

We have seen how to obtain the self-interaction correction to the probabil- 
ity of single particle emission. The natural next step would be to obtain 
similar results for multi-particle processes, involving some combination of in- 
going and outgoing particles. The ultimate goal, of course, it to construct 
a complete S-matrix relating arbitrary in and out states. This brings to the 
fore what is usually regarded as the central conceptual puzzle of black hole 
physics: does such an S-matrix exist which unitarily relates states described 
on J + and or is there "information loss" is the sense that pure states on 
J~ can evolve into mixed states on J + l While the latter possibility clearly 
violates a tenet of quantum physicists, there is at present no satisfactory 
proposal as to how the former possibility can be realized. 

If information is preserved in gravitational collapse, and an S-matrix ex- 
ists, it will require the existence of intricate correlations between particles on 
J + which encode the details of the state on J~ . An understanding of how 
this situation might arise is currently precluded by the fact that essentially 
nothing is known about how to compute any correlations on J + , much less 
those which would preserve all information. 

In free field theory, the state on J + is described by an exactly thermal 
density matrix, and no one knows how to go beyond the free field approx- 
imation, except in the case of the single particle processes we have been 
discussing. However, we can envision calculating the correlations between 
two emitted particles by an extension of our previous methods, simply in- 
cluding two shells instead of one. Presumably, the correlations would be quite 
complicated in the case of short time separation between the particles, but 
upon going to to large time separations one would see the later particle being 
emitted with a probability corresponding to a hole of mass M — u, where cu 
is the energy of the earlier particle. In addition, there is the possibility for 
the phases of the particles to be correlated even for large time separation, 
owing to the fact that the two shells were closely spaced near the horizon at 
early times. 

There is no obstacle to writing fown the action for two shells, 

S = Jdt [pir'i + p 2 r 2 - H(r 1 ,r 2 ,pi,p 2 )]. (4.105) 

H is again given by the ADM mass, which can be expressed in terms of the 
shell variables by solving the constraints. The action S^k^t, r±, r 2 ) can then 
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be computed for the initial condition 

£fcifc 2 (0, ri, r 2 ) = hri + k 2 r 2 (4.106) 

by integrating the Hamilton- Jacobi equation in precisely the same manner as 
for the one particle case. Presumably, all correlations between the particles 
are encoded in this action; unfortunately, we do not know the code, for 
reasons that we now discuss. 

In the single particle case it was possible to make progress because we 
know how to pass from the first quantized description in terms of Sk(t, r), to 
a field description; namely by writing 

0(t, r) = J2 [a k e lSk{t ' r) + ale~ iSk{t ' r) ]. (4.107) 

k 

This then put the full power of Bogoliubov transformations at our disposal, 
which enabled us to determine the emission probabilities. By constrast, in 
the multi-particle case the relation between the first quantized description 
and the field description is unclear — we are not aware of any extension of 
fl4.107| ) which takes into account the two particle solutions. This represents 



a major obstacle, since a proper interpretation of the theory requires that we 
have a field description. On the other hand, it seems most likely that this 
problem is merely a technical one, as there is no reason to believe that the 
theory is ill defined or inconsistent. 



Chapter 5 

Black Hole Entropy 



Almost all researchers agree that a black hole has an entropy equal to one 
quarter the area of its event horizon, even though there is no consensus as 
to what the entropy represents physically. The most appealing possibility 
is that the entropy counts the number of black hole microstates, i.e. the 
black hole Hamiltonian has e s( - M ^ eigenvalues between and M. Since it 
may seem rather remarkable that this could be deduced from the study of 
free field theory on a classical background geometry, we will attempt put the 
argument in perspective by recalling the corresponding derivation in ordinary 
thermodynamics. Suppose we wish to determine the entropy of some body 
of matter. To proceed, we would first perform the experiment of placing the 
body in contact with a heat bath of temperature T, waiting for equilibrium 
to occur, and then measuring E(T), the energy of the body. Then, assuming 
that equilibrium occurs when the total number of states of the system is 
maximized, we have 

dS hody (E) _ dSbath 

dE ~~dE~- (5 - lj 
Finally, since dS hath /dE = 1/T, we obtain 

= fiw (5 - 2) 

Clearly, the derivation relies on two features: a) our ability to measure 
E(T), or T(E), b) the assumption that in equilibrium the number of states 
is maximized. Now let is return to the black hole case. Obviously, no one 
has performed an experiment to see whether a black hole of mass M is in 
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equilibrium with a heat bath of temperature l/8irM. The point is that it 
appears that any nonsingular state of the black hole satisfies this property. 
We saw that by requiring regularity of the stress energy tensor at the horizon 
we were inevitably led to consider states which radiate at the Hawking tem- 
perature. The status of point (b) is much less clear. Just because the hole is 
in equilibrium with the heat bath we cannot conclude that the black hole is 
making frequent transitions among its microstates, and that the probability 
distribution of these states is Boltzmann. A rather trivial example of this 
behaviour is provided by a perfectly reflecting mirror; it can certainly be 
in equilibrium with thermal radiation without satisfying these other condi- 
tions. In fact, we shall se shortly that it does not seem to be possible to 
assign a Boltzmann probability distribution to the microstates of the hole. 
Nevertheless, let us assume for now that assumption (b) is justified. Then, 

r M dM A 

The belief that the entropy does count microstates is bolstered by the 
existence of a completely unrelated derivation due to Gibbons and Hawking 
[ |33|j . We will not go through the details but simply sketch the idea. Their 
approach relies on the fact that the partition function for a system can be 
expressed as a path integral over configurations periodic in imaginary time: 

Z = Tre-f )H = ( T>g^e- S . (5.4) 

J periodic 

Here, the path integral is over Euclidean metrics with periodicity Ate = i(3. 
One then calculates the path integral in saddle point approximation and 
notes that only the analytic continuation of the black hole geometry with 
f3 = 8nM contributes, because other geometries have a conical singularity 
at the horizon. After calculating the action and subtracting the flat space 
contribution, they obtain: 

Z{p) = e A/A e- pM (5.5) 

leading to the identification S = A/ '4. 

This derivation also has serious problems, not least of which is the fact 
that the Euclidean path integral for gravity may not even exist due to the well 
known conformal instability. Nevertheless, it has the virtue that the answer 
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agrees with the previous result. In any event, it seems overwhelmingly likely 
that the entropy has some sort of deep significance, and the most natural 
interpretation is that it counts microstates. 



5.1 State Counting 



If we are inclined to believe that a black hole has the enormous number of 
states e s , it behooves us to explain how these states can be understood in 
terms of the underlying Hamiltonian. It would be most satisfying if quantiza- 
tion of this Hamiltonian revealed a discrete set of states whose number could 
be counted to yield the entropy. In this section we will describe some efforts 



along these lines, and the divergence problems |34| which ensue. In keeping 
with the spirit of this thesis, the discussion will be confined to spherically 
symmetric configurations. While we have no particular reason to believe that 
all of the states can be accounted for by considering only the s-wave, the sorts 
of problems that arise seem to afflict the higher partial waves in the same 
way. Pure gravity in the s-wave has no dynamical degrees of freedom, and 
so to have a nontrivial theory to quantize we include a massless scalar field. 
It is a little unsettling that we are forced to include to matter in order to 
calculate, since the expression for black hole entropy has no explicit depen- 
dence on the matter content of the world. However, it is possible that there 
is matter dependence, but it can be absorbed into the value of Newton's 



constant which appears in the entropy formula [36. 37 1 



As a first approach, we can try to proceed in the same way as if we were 
quantizing a soliton in flat space — by quantizing the quadratic fluctuations 
of the field about the background solution. Considering only the quadratic 
fluctuations amounts to doing free field theory. For simplicity, we will work in 
Schwarzschild coordinates and use modes of definite energy, e - luJ ( t - r *)^ where 
r* = r + 2mlog(r/2M — 1), and write: 

= EK e " Mt_r * ) +4e M *- r * ) ]. (5.6) 

We consider only the outgoing modes since they alone lead to the divergence 
problems. In order to count modes we must make the frequencies discrete in 
some way. The easiest thing to do is to impose periodic boundary conditions 
at r = 2M + e and r = L. The e is included because the modes become 
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singular at the horizon, and L is some arbitrarily chosen radius outside the 
horizon. We take L ^> 2M ^> e. This leads to the allowed frequencies: 

Un * M log (LA) ' (5 ' 7) 

This means that there are 

n = M "'° g(L/6) (5.8) 

7T 

single particle states for the black hole with energies between M and M + oo. 
The problem is that there is no physical reason to keep e finite, but the 
number of states diverges as e — > 0. The reason for this behavior is due to 
the fact that — > — oo as r — > 2M, so the modes oscillate an infinite number 
of times before they reach the horizon. 

It seems quite possible that the infinite number of oscillations is simply 
a result of using free field theory, and that once the proper self-interaction 
corrections are included a finite result will be obtained. To illustrate this, 
let us recall the expression for the canonical momentum of a self-interacting 
particle, 



Pc 



V2Mr - J2M + r - r log 



For an energy eigenstate M + = M + u, and the corresponding mode is 



(5.9) 



il>(r,t) = e i f rp ^ r ' )dr '- iu ' t . (5.10) 

Thus the number of oscillations is finite or infinite depending on whether 
J 2 ( M+u; ) p c (r')dr' is finite or infinite. The singular part of p c as r — > 2M is: 

p c ^ -2(M + cu)log[r -2(M + cu)] (5.11) 

which leads to, 



/ 



2{M+uj) 

Pc (r')dr' = finite. (5.12) 



The (incorrect) free field theory result is recovered by expanding in uj. Then 
p c has the singular part 

AMuj 

p c = 5.13 

F r-2M y ' 



CHAPTER 5. BLACK HOLE ENTROPY 



69 



and J 2M p c dr is again infinite. 

Thus by including self-interaction we seem to have solved the divergence 
problem. However, there is a major caveat. The mode solutions ( 5.1U| ) were 
obtained using the WKB approximation. As we discussed earlier, the WKB 
approximation for a solution of the form e is only valid provided 



dS_ 

dr 



> 



d 2 S 



dr" 1 



d 3 S 



dr 3 



(5.14) 



In other words, S should be rapidly oscillating but the rate of change of 



oscillation must not be too large. Applying this condition to ( |5.10|) gives 



\Pc(r)\ > 



dp c {r) 



dr 



or: 



|2(M + u) log [r - 2(M + u)]\ > 



2(M + uj) 



(5.15) 



(5.16) 



r-2(M + tu) 

This condition is not satisfied as r -> 2(M + to). Therefore, we cannot 
trust the behavior of fl5.1U| ) near the horizon, including the conclusion that 
it oscillates a finite number of times. The most convincing resolution would 
be to go beyond the WKB approximation and obtain the correct result for 
if). But, as we have discussed previously, this requires the resolution of factor 
ordering problems which we do not know how to solve at the present time. 
The most we can say is that the preceding analysis suggests quite strongly 
that the divergences are connected with an improper treatment of the self- 
interaction of the modes near the horizon. 

The breakdown of the WKB approximation occurs when we insist on 
using modes of definite energy. Earlier, we saw how to derive a complete set 
of modes which were nonsingular at the horizon and are accurately described 
by the WKB approximation. These were: 



u k (t,r) 



JS k {t,r) 



(5.17) 



where: 

S k (t, r) = -(2M 2 - r 2 /2) log [1 + e W™'-*W)\ (5.18) 

with M' given by (|4.40| ). The coordinates are those of ( p.3|) . Since these 
modes do not have definite energy, there is no straightforward way to count 



CHAPTER 5. BLACK HOLE ENTROPY 



70 



states using the microcanonical ensemble. However, we can imagine putting 
the black hole in contact with a distant heat bath at the Hawking tempera- 
ture. If the black hole behaved as an ordinary thermodynamic body we could 
proceed by computing the partition function, and from that extract the en- 
tropy by standard thermodynamic formulas. Since the partition function is 
a trace, 

Z = Tr(e~ f>H ), (5.19) 

it is independent of which basis we choose to describe the states, and so the 
modes ( 5.17|) are as good as any other. Our working assumption is that we 
can describe all the states as fluctuations about a fixed black hole of mass 
M. Since in the s-wave there are no purely gravitational fluctuations, the 
Fock space built on the modes (|5.17|) should provide a complete description 
of these states. Therefore we can write 

Z =< 0|e _/3H |0 > + (dki< hle'^lh > + Jdh dk 2 < k^e'^lk^ > + ... 

(5.20) 

H is the total energy as measured at infinity, so e~^ H generates translations 
in imaginary time: 

e pH <j>(t)e- pH = <t>{t-i(3). (5.21) 
This feature makes Z easy to calculate. The first term gives simply 

< 0|e _/3H |0 >= 1 (5.22) 

by definition of the vacuum. The next term is more interesting, 

1 /"OO 1 roc 

< k\e~ pH \k >= — / dru* k (0,r)u k (-i/3,r) = — / dr e~ lSk{0 ' r) e iSk{ - if3 ' r) . 

2% J-oo 271 J-oo 

(5.23) 

Let us concentrate on the behavior for large k. For fc>Mwe have 

S k (-i/3,r) w 2M(r - 2M)(k/2M' + i(3 /AM') ; M'«M + ^(r-2M). 

(5.24) 

This is strictly valid only near the horizon, but that is the only region in which 
we need the solutions since we can always construct wavepackets localized 
near the horizon. The point we now wish to stress is that < k\e~ l3H \k > goes 
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to a (nonzero) constant independent of k as k — > oo. The precise value of 
the constant depends on the form of the wavepackets, but at any rate 



so that Z cannot be defined. 

There is a simple physical explanation for this behavior. A nonsingular 
mode has positive energy density for points outside the horizon, and neg- 
ative energy density for points inside. Now imagine being at fixed radius, 
r, initially outside the horizon, and letting k increase. As k increases, the 
total mass inside radius r increases. However, the mass can never increase 
past M = 2/r, since if it did one would be inside the horizon, but the modes 
have negative energy here and so cannot lead to an increase in mass. So as k 
goes to infinity, the effect on the geometry simply goes to a constant, which 
explains why the matrix element of e~@ H also goes to a constant. 

The lesson to be learned from all of this is presumably that a black hole 
cannot be treated as an ordinary thermodynamic body in the sense that 
its states are distributed according to a Boltzmann distribution when in 
equilibrium with a heat bath. There are simply too many low energy states 
localized near the horizon for this distribution to make sense. Let us point 
out, though, that this conclusion is based on the assumption that the states 
are correctly described by a local quantum field theory. This assumption may 
be incorrect, and the divergences may disappear when the correct theory at 
short distances, such as string theory, is taken into account. That string 
theory plays a crucial role in black hole physics is argued in, for example, 




(5.25) 



Ref. |3§. 



Chapter 6 

Black Holes and Quantum 
Tunneling 



The radiation of particles from matter evolving along a classical trajectory 
has been heavily studied in recent years. Less well studied is the radiation 
accompanying quantum tunneling from one classically allowed trajectory to 
another. The following question is of interest: if a matter system impinges 
upon a potential barrier with a radiation field in a certain state, what is the 
state of the field given that the matter is subsequently observed to be on the 
other side of the barrier? A method to answer this question in the context of 



false vacuum decay in flat space was developed by Rubakov |38j and has been 



generalized to include gravity as well as topology changing processes |39|, [40| 



The spectrum of radiation is found by solving an imaginary time Schrodinger 
equation, the occurrence of which leads to novel features. Instead of solving 
field equations in real time, one is naturally led to consider propagation on 
the Euclidean solution interpolating between the two classical trajectories. 
As phase factors in real time are converted into exponential damping factors 
in imaginary time, the resulting particle creation can be distinctly different 
and is accompanied by the systematic supression of excitations present before 
tunneling. 

Given this situation, it is natural to ask how the radiation from black holes 
might be affected by the presence of tunneling. If we consider a distribution 
of matter, initally outside its Schwarzschild radius, which tunnels through 
a potential barrier to form a black hole, the conventional calculation of 
the radiation does not apply. On the other hand, it would be shocking if 
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the same answer was not obtained for the radiation at late times, as this 
is thought to depend only on the hole's late time geometry and not on its 
history at early times. Here we compute the radiation for this process and 
show that while the Euclidean time evolution has an effect at early times, it 
has none at late times so that the standard result is in fact obtained. 



In order illustrate the technique of Ref. [38] in a simpler setting, we first 



study the effect of tunneling on another well known radiating system — the 



moving mirror pBfl . We show in Sect. (2) how an imaginary time Schrodinger 
equation emerges from a Born-Oppenheimer approximation, and use this 
result to calculate the shift in the spectrum of radiated particles as a result 
of the tunneling. It is shown that the initial spectrum is shifted to favor 
low energy excitations, as is understood by realizing that the probability to 
tunnel is increased if energy is transferred from the radiation to the mirror. 

In Sect. (3) this approach is extended to include gravity in asymptotically 
flat space. A WKB approximation to the Wheeler-DeWitt equation, as con- 
sidered in f4l| , 42] , is used to obtain an imaginary time Schrodinger equation 



which can then be solved as before. In Sect. (4) we use this result to examine 
the radiation from a black hole which is formed by tunneling. In particular, 
we consider the tunneling of a false vacuum bubble, a system extensively 



studied in Refs. [|3[] - [|7|. This example involves a complication due to the 
peculiar structure that arises; Refs. [(48], [Tj| show that the sequence of three- 
geometries encountered during tunneling can not be stacked together to form 
a manifold. Employing a slight modification of the standard approach, we 
show how the behaviour of fields on the Euclidean Schwarzschild manifold 
protects the late time radiation from being affected by tunneling. An intu- 
itive reason for this is that the bubble's tunneling probability is unchanged 
by the presence of Hawking radiation, which involves the creation of pairs of 
particles with zero total energy. 



6.1 Tunneling Mirror 

Consider a mirror moving in a one dimensional potential in the presence of 
a massless scalar field. The Schrodinger equation for this system is 

d 

[H m + H^[<j),x m ;t} =i— tf[0,x m ;t] (6.1) 
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where 



and 



Ha, 



H„ 

1 

' 2 



1 d 2 
2m dx£ 



dx 



5(j){x) 



+ V{x m ) 



+ [-r- 



dx 



(6.2) 



(6.3) 



Note that \l/ is a function of the mirror coordinate x m , and a functional of 
the field configuration <f>{x). The mirror boundary condition is imposed by 
demanding that the field vanish at x m , 



V[<p,x m ;t]=0 if 



(6.4) 



The system is solved by assuming that the backreaction of the field on the 
mirror is a small perturbation of the mirror's motion, and that the mass 
and momenta of the mirror are large enough that it can be described by 
a well localized wave packet. In this domain the system admits a Born- 
Oppenheimer approximation, which amounts to an expansion in 1/m. In 
particular, we seek a solution to the time independent Schrodinger equation 



[H m + [0, X m ] = E^[4>, X r 



(6.5) 



valid to zeroth order in 1/m. Following Refs. 0, 42] the Born-Oppenheimer 
approximation is implemented by writing \1/ in the form 



X m \ = 1pw(Xm) e iS{Xm) X[<t>, ^ 



(6.6) 



where ipvv is a slowly varying function to be identified with the Van Vleck 
determinant. To lowest order in 1/m, ( |6.5|) reduces to the Hamilton- Jacobi 
equation. 




V(x r 



since dS/dx m , V(x m ) and E are all of order m. 
To zeroth order: 



E 



(6.7) 



i d 2 S 
2m dx m 

dS 



— ^Pvv-, — ^— 



ipvvX[<P,x 

d 



dS dip 



vv 



rn doc doc y. 



(6.8) 



x m \ + ip vv H^ x[<f>, x Tl 



0. 
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ipvv is chosen so that the first two terms cancel, leaving 



H,f>x[<i>,x m ] = d x[0,a; ro ]. (6.9) 



This can be put in a familiar form by defining the time variable r(x m ). In 
a classically allowed region, where E — V(x m ) > and dS/dx m is real, r is 
defined by 

- — = — allowed regions (6.10) 
whereas in a classically forbidden region with dS/dx m imaginary, 

— — — = i — forbidden regions. (6-H) 

The resulting zeroth order equations for are: 

d 

H<f> XWi T ] = i~Q- r ] allowed regions (6-12) 
d 

- Hcj, x[(f>, t e ) = -k — x[<t>, t e] forbidden regions. (6.13) 

These are the fundamental equations governing the evolution of the scalar 
field in the presence of the mirror. In the allowed regions we have recovered 
the time-dependent Schrodinger equation with the postion of the mirror play- 
ing the role of a clock, whereas in the forbidden regions we have obtained 
a diffusion equation, which we shall refer to as the Euclidean Schrodinger 
equation, with the Euclidean time te measuring the position of the mirror 
in the potential barrier. 

Now, choose the potential to be of the form illustrated in Fig. |6.1| and 
let the mirror come from right to left. In the allowed region to the right of 
x l m the state x[(p, T ] obeys the normal Schrodinger equation, and so standard 
methods can be used to find t 1 }- Between x l m and x^ the mirror is in a 
forbidden region, so the state evolves according to 



— / dx 

2 Jx m (T E ) 



d 2 



5<j){x) 




Atb] = ^-x[0,Te] (6.14) 
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V(x m ) 




Figure 6.1: A generic mirror potential. The turning points for energy E are 
indicated. 



with xI^iTje] = x[4>i T% ]- We wish to solve this equation in order to find the 
state at the final turning point x^. It is useful to transform the mirror to 
rest by defining the coordinate 



y{x,T E ) = x - x m {r E ) 
in terms of which the Euclidean Schrodinger equation is 



-5JT* 



5 2 dx m deb 5 

+ 2—^--—— + 



s<l>(y) 



dr E dy5(f)(y) \dy 



x[<t>,T E \ 



_d_ 

drf 



or 



The solution is 



x[<P,r E ] = Texp 



-H^r E )x[^r E } 



_d_ 



x[(f>,r E } 



(6.15) 



X[0, T E \ 

(6.16) 
(6.17) 



I" Hf(r E )dT E 

Jt'„ 



x[<P,t e ] = U e (t e ,t e )x[<P,t e ]. (6.18) 



Here T represents time ordering with respect to r E . The crucial point is that 
the Euclidean time evolution operator, U E , is non-unitary. This is natural 
since we know that wavefunctions decay exponentially during tunneling. If 
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Ue was unitary, the easiest way to calculate it would be to transform to the 
Heisenberg picture, solve the field equations mode by mode, and compute 
Bogolubov coefficients. However, as emphasized in Ref. [38] the non-unitarity 
of Ue implies that the Schrodinger and Heisenberg pictures are inequivalent, 
making the standard method inapplicable. Instead, one can use the method 
developed in Ref. |3pJ which closely resembles the standard one but is more 
general. We first describe the state right before tunneling. For convenience, 
set x % m = t % = t e = 0. Let £ w (x,t) be a complete set of positive norm 
solutions to the Klein-Gordon equation which vanish vanish at the mirror: 



d 2 



+ 



o 2 



dr 2 dx 2 



€w{x,t) = 



dx 







d 

OT 



(6.19) 



(6.20) 
(6.21) 



The set of allowed frequencies u is taken to be discrete, and J2u represents 
summation over this set. The field operators can then be expanded in terms 
of these modes: 



with 



TT<p(x,r) = 



UJ 



(6.22) 
(6.23) 



Now define Euclidean fields <fi E (y, te), n E (y, t e) which agree with <f>(x, r), 



it^y, r) at r = r = 0, but evolve according to 

s (y, t e ) = Ue 1 (t e , 0) 4> E {y, 0) U e (t e , 0) 
n E (y, t e ) = U e \t e , 0) 7t E (y, 0) U e (t e , 0). 

We will calculate Ue(te,0) by first finding cf> E (y,TE), ^s(Vi t e) 
equations for these operators are 



8tf 



dr E dy 



(6.24) 
(6.25) 
The field 

(6.26) 
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.d 2 4> E dx m dit* 
-i - „ + 



dr E dy 



So 



and 



2aE 



dri 



+ 



drr, 




dr E dydr E dr£ dy 



Equation ( |6.29| ) can be obtained by varying the action 

S = 2 / dydTE 9 z d n 
with the Euclidean metric 



ds E = g E dx^dx v = dr^ + 2 — — dx dr E + dx 

dr E 



(6.27) 
(6.28) 

0. (6.29) 

(6.30) 
(6.31) 



<f) E , it E can be expanded in terms of modes f w which satisfy the Euclidean 
Klein-Gordon equation (|6.29|) and which vanish at y = 0, 

t E (y,T E ) = J2bu J Uy,T E ) (6.32) 



Trf (y, tb) = iJ2 ^ 



d , , dx m d 
■fw(y,T E ) - 



fu>(y,T E ) 



(6.33) 



dr E ' dr E dy 

As the Euclidean Klein-Gordon equation is elliptic, one cannot in general 
impose Cauchy boundary conditions at t e = on f u . The resulting solu- 
tions would not satisfy the mirror boundary condition. With the appropriate 
boundary conditions, either Dirichlet or Neumann, imposed at r E = and 
t e = t e , a detailed calculation is, of course, required to find f u for a generic 
mirror trajectory. We shall take the solutions as given and only use their 
specific forms in a region far from the mirror, where they are simple. 

Now, using the condition that the two sets of operators </>, and (fi E , n E 

are equal at r = t e = 0, and taking inner products, the operators can be 
expressed as a linear combination of a^, d^: 



'uiuj'O'u,' 



(6.34) 
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Then using 

t E (y, 4) = U e \t e , 0) My, 0) U E (r f E , 0) = U E \r f E , 0) 4>(y, 0) U E (r E , 0) 

(6.35) 

and the analogous expression for n^, the following equations for U E are 
obtained: 
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E E [awou* + Putf^ fu(y, t~e) 



E [Ue\4, 0) a w C/ B (r|, 0) £,(y, 0) + ^(rj, 0) a* ^(r|, 0) £(y, o/ 

a; 

(6.36) 



and 



E 



d 

i E E [owouy + « — / w (y, t|) 

^(rl, 0) a w ^(r|, 0) —&(!/, 0) + C/^(r|, 0) a+ C/ £ (r|, 0) — £(y, 0) 

(6.37) 

(6.38) 
(6.39) 



Again taking inner products, this leads to relations of the form 



Uz\Ti,0)a u U E (T f E ,0) = Y,[ 



U e 1 (t e , 0) al U E (r E , 0) = E [wuxSGu 1 + 



Then it can be shown that [38 
Ue{t e ,0) = const, x : exp^^ 



where the matrices D, F, and G are defined by 



D = vz- 1 ; F = 
The state after tunneling is then determined, 

X (4)) = u E (4) | x (o)> 



1 ; G = -Z- X W. 



(6.40) 
(6.41) 

(6.42) 



and is expressed in terms of occupation numbers with respect to the modes 
£u>(y> 0), where now y = x — x^ m . All of the information about the final state 
is contained in the matrices D, F, and G, which are in turn given in terms 
of inner products between the modes f u and 

As a simple application of these formulae we will calculate the shift in the 
spectrum of outgoing particles which are far from the mirror at the time of 
tunneling. It is assumed that the mirror was initially at rest and the field 
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in its ground state. The mirror subsequently accelerates in the potential 
V(x m ) until it reaches the classical turning point x l m . It is well known that 
as a result of the mirror's acceleration, a flux of outgoing particles is created 
whose spectrum is calculable by standard methods p8| . Outgoing particles 
far from the mirror are wavepackets composed of superpositions of plane 
waves, 

£,(x,r) = J_e-M™0 ( 6 .43) 

2y/UJ 

The spectrum of outgoing particles located at x = x 3> uj^ 1 at r = is 
written as 

E^(K})IK» (6-44) 

where {n^} is a set of occupation numbers and ({n^}) is the amplitude 
for the set to occur. 

Far from the mirror, the modes f u are easy to calculate since the mirror 
boundary condition is irrelevant. They are of two types, 

t~ _ \ ^-U}T E +iulX _ \ e -U)T E +iu(y+Xm(T E )) 



f + = _ e u)T E +iujx _ _ UT E +iu){y+x m (T E )) /g ^rt 

Ju 2^J ~ 2^J 1 ' ; 

Then 0, n and <fi E , are equal at r = r E = if 

K = Ow ; it = ai. (6.46) 



Equation ( |6.38| ) gives: 

U E \r E , 0) U E (4),0) = e — 1+-^ a u 
Ue'H, 0) at U E (r E , 0) = e ^+-^ at (6.47) 

leading to 

D = G = ; = (e-»4-*«>>L _ A ^ ( 6 . 48 ) 



and 



U e (t e ,0) = const, x : exp V] 



-UJT^-iuJxb _ ^ 



a},a„ 
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= const, x : expj^ e lu3Xm - 1 ala w :: exp^] e ut e — 1 c^Ja^ : (6.49) 

The first factor is a translation operator which expresses the state in terms 
of the modes £ u (x,0) instead of ^(x + x-^,0), and the second factor acts 
on a state |{^}) to give e~ E ^ nui '^ T B where E ({n u }) = Y. n uj UJ is the 

energy of the state. Therefore, the state after tunneling is 

const, x Y, e- £({n " }) ^ S s ({n u }) |K» . (6.50) 

{nui} 

The result of the tunneling is simply to shift the spectrum from S s to 

It is not difficult to understand this result. Since the total energy is 
fixed , the state before tunneling is given by a superpostion of the various 
ways of distributing the energy between the mirror and the radiation. As 
the mirror's probability to tunnel depends exponentially on its energy, we 
expect an inverse exponential correlation between tunneling and energy in 
radiation. Thus an observer measuring the spectrum of radiation, conditional 
on the mirror tunneling, finds the result ( |6.50| ). Far from the mirror the shift 
in the spectrum depends only on Tg, the amount of Euclidean time spent 
during tunneling. This is because the tunneling amplitude in the WKB 
approximation is e~ s , and the derivative of S with respect to energy is just 
the Euclidean time. 

If we were to identify the Euclidean time with an inverse temperature, 
the shift would become a Boltzmann factor. This makes it easy to gener- 
ate thermal distributions of radiation. Specifically, if the distribution before 
tunneling was a constant, then after tunneling tracing over the states of the 
mirror would yield a thermal density matrix for the radiation. A number of 
authors have been led by this fact to seek a connection between the thermal 
radiation that arises in the contexts of cosmology and black holes and an 



occurrence of tunneling [f40| [50], |5T|. Such a connection relies upon assump- 
tions about what is on the other side of the barrier and what the spectrum of 
radiation is there. In this work we only consider situations where there is a 
well defined classical trajectory on either side of the barrier; we are interested 
in the case in which there is collapsing matter on side of the barrier and a 
black hole on the other. The treatment of this process requires an extension 
of the previous method to include gravity. 
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6.2 Application to Gravity 

In this section we make a WKB approximation to gravity in a manner which 
directly parallels that for the moving mirror. The action for gravity plus 
matter takes the form 

2 

S = —r- / d^xJ^g (71 — 2A) + Sm + boundary terms 

167T J 

= J d A x (tx^ + nijW - N l l-L t - N{hC) + boundary terms. (6.51) 
The Wheeler-DeWitt equation resulting from this action is 



167r 5 5 m 2 



1 /o \ 



= (6.52) 



m 2 Shij 5h k i 167r 
Proceeding as before, we seek a semiclassical solution of the form 

[hij, <Pi\ = ^wihij] e im ? s ^ x [fa hij] . (6.53) 
At first order the Einstein-Hamilton- Jacobi equation is obtained: 

^ G ««^lr-^( 3K - 2A ) = - < M4 > 

m 2 ShijShkL 16tt V > 
Zeroth order yields 

16tt SS 5x t ^ ( v 
2 1 G ijkijr-j7- + W tM x = 6.55 

provided ipvv satisfies 

5 2 S . _ 55 5^ 



G vfci 7r k ^ + Gwiirifir- = °- ( 6 - 56 ) 



'ShijShij 5hij 5h kl 

The momentum constraints at first order are 

/ 55 



(6.57) 
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and at zeroth order are 



2? 



$X 



Sh 



(6.58) 



v / \j 



Equations ( |6.55| ) and (|6.58| ) describe how the matter wave function evolves 
as the spatial geometry changes. Quantum field theory in curved space can 
be recovered by writing x's dependence on hij in terms of a time functional 
t[x; h^}, and by reintroducing a lapse N T and shift N iy demanding that they 
obey 



ijkl 



SS_ 

Shu 



m' 



16ttN t 



dy 



Sh 



ki 



Then 



16tt 



— i- 



N T G ijkl -^-^ + 2iN i 
Shij Sh k i 



Sr[y; h ab ] 



Sh 



v ' li, 



Shij Sx 



St Sh 



(6.59) 



(6.60) 



it 



so that the equation for \ becomes 



drx 



N T H tM + N.H 



xl 



i-^xifih r ] 



(6.61) 



The condition ( |6.59| ) agrees with the classical relation between and h^, 
demonstrating that r[x\ hij] is the classical time and that ( |6.61| ) is the Schrodinger 
picture version of quantum field theory in curved space. 

As with the mirror example, r becomes imaginary during tunneling so 
we define a Euclidean time te along with a Euclidean lapse N TE = iN T , in 
terms of which \ obeys 



- / dx[N^H tM + iNiH^xtti, r E ] 



d 
dr E 



x[&>Te] 



For a massless scalar field with action 



we have 



n 



tM 



hr*Tti + Ksh %3 di6di 



(6.62) 

(6.63) 
(6.64) 
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n 



(6.65) 



To evolve \ through the tunneling region one is required to calculate the 
Euclidean time evolution operator 



Ue(4,t e ) 



Texp 



H^dr E 



(6.66) 



with 



Hi 



I dm I om 



(6.67) 



As before, one proceeds by defining Euclidean fields obeying ( |6.24| , |BT25l) . In 
the present case the resulting field equations are: 



7rf = % 



. hh (da 



N TE \dr E 







N l di 



(6.68) 



with 



dsl = g^ u dx^dx u = (N TE dr E ) 2 + hij (dx* + NidT E ) (dx j + N j dr E ) . (6.69) 

The evolution operator, and therefore the state after tunneling, is determined 
by solving the field equations mode by mode, and repeating the steps leading 
from (0|) to 



6.3 Black Hole Radiance in the Presence of 
Tunneling 

We can now apply this method to determine how the radiation from a black 
hole is affected by tunneling. It is well known that a black hole formed 
classically from collapsing matter radiates in a complicated manner at early 
times due to the time dependent geometry, but at late times will inevitably 
radiate as a black body at the Hawking temperature. Is this scenario altered 
if the black hole is formed while tunneling? We shall show that it is not. The 
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form of the late time radiation is insensitive to the hole's unconventional 
history in a way that is consistent with the intuitive picture of Hawking 
radiation being caused by pair production near the horizon. 

We consider the behaviour of a scalar field on the background of a false 
vacuum bubble which tunnels leading to the formation of a black hole. The 
action for a false vacuum bubble in the thin wall approximation is 

S = — f d^x^U - — / d^xV^F -—I d?A (6.70) 

167T J 07T Jbubble 4.7T J wall 

where Aj is the cosmological constant of the false vacuum, and \x is the 
energy density of the bubble wall. The classical solutions for this action have 
been derived in Refs. [fi3"|-[47|. In what follows we refer to the treatment 



of Ref. p7| |. The spherically symmetric solutions are characterized by three 
parameters: A/, /i, and the total mass M. In addition, for given Aj and \i 
there is a critical mass M cr below which there are two solutions: type (a), 
where the bubble emerges from a singularity with zero radius, subsequently 
expands to a maximum radius, and then recollapses; type (b), where the 
bubble initially collapses from infinite radius, reaches a minimum radius, and 
then reexpands. Using the results of Refs. j|8], [Tj| , we focus on an expanding 
solution of type (a) which tunnels to an expanding solution of type (b). We 
confine our interest to the region outside the bubble where the metric, written 
in terms of Schwarzschild time t and r* = r + 2Mln(r/2M — 1), is 

As t and r* cover only part of the complete manifold, we introduce Kruskal- 
Szekeres coordinates, 

ds 2 = (-dT 2 + dX 2 ) + r 2 dVL 2 . (6.72) 

r 

The two sets of coordinates are related by 



/ r 



A e r/2M =X 2 



\2M 

4Mtanh" 1 (T/X) if \T/X\ < 1 
4Mtanh~ 1 (X/T) if \T/X\ > 1 



(6.73) 
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Figure 6.2: The type (a) and (b) solutions. The heavy lines represent the 
bubble trajectory, and the dashed lines are the initial and final surfaces of 
the tunneling solution. In these figures, only the regions to the right of the 
trajectory are of interest, as they are outside of the bubble. 



Using these cordinates the type (a) and (b) solutions of interest are depicted 
in Fig. B72. 



The tunneling amplitude for this process has been computed by two dif- 
ferent methods. In Ref. [[TJ| the solution to the Wheeler-DeWitt equation is 
found in the WKB approximation by solving the Einstein-Hamilton- Jacobi 
equation ( 6.54Q . Since the solution behaves as e~ s , and the tunneling am- 



plitude is given by the ratio of the wavefunction evaluated at the initial and 
final geometries, the tunneling amplitude is 



exp (S[h 



initial] _ ^flnal]) ( g ?4) 



No difficulties arise in this approach; the calculation of tunneling amplitude 
proceeds in a straightforward fashion. 



In Ref. [|8[ the calculation is performed using the functional integral. In 
this formalism one looks for a manifold which interpolates between the initial 
and final surfaces and which is a solution to the Euclidean Einstein equations. 
The tunneling amplitude is e~ s , where S is the action of the solution. It is 
found, however, that solving the field equations leads to a sequence of three 
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geometries which do not form a manifold. To see this, first note that the 
geometry outside the bubble is Euclidean Schwarzschild space, obtained by 
t -> it E , T -> iT E , 



ds\ 



with 



2M\/ l2 , 2 n - a 32M 3 e- r / 2M / - , ~x , - 
I ] ( dt\ + dr 2 ) + r 2 dn 2 = (dT 2 + dX 2 ) + r 2 dtt 2 



r 



r 



2M 



lj e r/2rn = X 2 + T% ■ t E = AMtMT l (T E /X). 



(6.75) 



(6.76) 



It remains to describe the motion of the bubble wall. Solving the equations 



of motion leads to the trajectory in Fig. 3.3 



. final surface 



4- 




X 



Figure 6.3: Bubble trajectory in Euclidean Schwarzschild space. 



It is seen that the bubble wall crosses the initial surface during the course 
of its motion, creating a situation in which it is impossible to identify a 
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region which is swept out by the evolving hypersurface. Some regions of the 
manifold are crossed twice by the hypersurface, some once, and some not 
at all. The authors of Ref. call this object a pseudomanifold and give 



a prescription to calculate its action by assigning covering numbers to the 
various regions, but this is not needed for what follows. 

With these results in hand, the technique of Sect. (3) can be used to 
calculate the state of the scalar field after tunneling. It was seen that once 
the solution of the Einstein-Hamilton- Jacobi equation is given, the field wave 
functional x is fully determined by equations (|6.55| ) and (|6.58|) . Since S[hij] 



is calculated in Ref. WB\, we have all that we need to find x- This would, 
however, require finding the solution to an unfamiliar functional differential 
equation. To cast it in in the form of the Schrodinger equation a lapse 
N TE , shift iV, and time te were reintroduced leading to the appearance of 
the Euclidean metric g®. In the present case there is no true interpolating 
Euclidean manifold, so that any choice of N TE and iVj which define a well 
behaved will lead to a bubble trajectory that is a multivalued function of 
time. Alternatively, a choice of time functional which gives a single valued 
bubble trajectory will necessarily lead to a Euclidean metric with vanishing 
determinant at some point. In either case, it is not clear that the resulting 
Schrodinger equation is well defined. This is apparent from Fig. |6.3| , where it 
can be seen that boundary conditions imposed on the initial surface and on 
the bubble wall may contradict each other. These difficulties arise as a result 
of trying to compute the final state of the field in one step, which requires 
a Euclidean manifold interpolating all the way from the initial surface to 
the final surface, and can be avoided by calculating the state on a series 
of intermediate hypersurfaces. In this approach, it does not matter that 
the bubble wall eventually crosses the initial surface since once the state is 
calculated at some intermediate point we can forget about what preceded it. 

For simplicity, we will consider only the s-wave component of the scalar 
field and frequencies high enough such that the geometrical optics approxi- 
mation is valid. This means that the field equation is taken to be 

ggdpdrf = 0. (6.77) 



The state of the field on the initial surface, t = T = 0, is most conveniently 
expressed in terms of the coordinates r* and t. We divide the modes into 



CHAPTER 6. BLACK HOLES AND QUANTUM TUNNELING 



90 



ingoing and outgoing, 

&(r.,t) = C u e-**» r J 

e: t (r^)=ae" iw(t - r,) (6.78) 
and write the field operator as 

fa; t) = E [«L n C + a^£*(n, t) + in - ont] . (6.79) 

C u are normalization constants whose values will not be important. We shall 
only consider the in modes as the treatment of the out modes is exactly the 
same. We also suppress the in superscript. 

In the first stage of the evolution the hypersurface is pivoted around 
= r\ by 180°, where is the position of the bubble wall on the initial 
surface. The solutions to the Euclidean field equations are most conveniently 
obtained by choosing Cauchy boundary conditions on the initial surface, 
(clearly a valid procedure in this case) 

fi(r.,0)=Ur.,0) ; ^/>*,0) = -^C(r*,0) 

/«7(r.,0)=£(r„0) ; — f~(n, 0) = -<^£(r., 0) (6.80) 

It is also easiest to use the X, T coordinates as they are well behaved every- 
where. Since the evolution of the hypersurface is simply a reflection about 
the point X = X b , a mode which has the form f(X, T E ) on the initial surface 
has the form f(—X + 2X b ,T) on the new surface. Using the relations 

r, = AMln^X 2 + T E ; t E = 4Mtan~ 1 (T B /X) (6.81) 

and that 

f»(r*,t E ) = C u e ±ut *- i » r - (6.82) 
near the initial surface, one sees that near the new surface, 

f±(X, T E ) = C„ exp ( j^j^ - 4iM W ln(-X + 2X 6 )) . (6.83) 

Since on the new surface, /+ = (/j)* and df+ /dt E = —(df~/dt E )*, the 
evolution operator U E is unitary. This means that the state on the new 
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surface has the same form as it did on the initial surface, but is now expressed 
in terms of the modes 

UX,T) = C u exp ( ~x + 2X» + 4tMujln (- X + 2X ^) " ( 6 ' 84 ) 

These modes can be approximated near T = as 

_ f C B e^.) if \X\»X b 

\ C u e -( 2iM ^/ xb )( T - x ) if|X|<X 6 1 J 

Now it is useful to express the state in terms of modes which are nonzero 
only inside or outside the horizon, 



D ^ e Mt-r t ) if X < 

if X > 




uxtl. < 6 ' 86) 

A fundamental result || in the derivation of black hole radiance is that 
the vacuum state with respect to modes which have a time dependence e~ lwT 
is the state 

const, x J2 e~ E ^/ 2T « \{n UJ }} < \{n (J }) > (6.87) 

{n u } 

with respect to the modes r]< and rj>. The sum runs over all sets of occupa- 
tion numbers, E = J2 n w^, arid Th = l/87rM is the Hawking temperature. 
Further, near the horizon, any deviation of \\) from the vacuum state can 
be ignored because of the arbitrarily large redshift as — > — oo. Far from 
the horizon £ w and rj^ agree so the form of the state is unchanged there. 

Now the hypersurface can be evolved the remainder of the way. If we 
restrict our attention to the region X < X b , then the motion of the hy- 
persurface is simply a translation, t E — — At#. This causes states with 
time dependence e luJt to be damped by a factor e _wA * B , and states with time 
dependence e~ luJt to be amplified by a factor e wA * B . Near the horizon, the 
state \x) consists of pairs of positive and negative frequency states according 
to ( |6.87| ). One member of the pair is damped but the other is amplified by a 
compensating amount so as to leave the state \x) unchanged. The final state 
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of the field can then be summarized as follows. Far from the hole, where 
there is no pairing, the initial state is damped: 

]T |K}> — const, x £ e~ E ^ AtE S({n w }) \{n w }) . (6.88) 

{riuj} {n^} 



Near the horizon the final state is given by (|6.87| ). This is true for both the in 



and out modes, so an observer stationed on either side of the horizon would 
observe a thermal distribution of both ingoing and outgoing particles. As 
time passes, all of the ingoing particles will eventually cross the horizon and 
be swallowed by the hole, whereas the outgoing particles will propagate out 
to infinity where they can be detected at arbitarily late times as a flux of 
thermal radiation at the Hawking temperature. 



6.4 Comments 

It was shown that the standard picture of black hole radiance is unchanged 
by tunneling. At late times, the hole radiates just as it would have had 
it been formed from a classical collapse. This makes sense if one thinks of 
Hawking radiation as pair production. The probability of tunneling is not 
affected by the creation of a pair, since the pair has zero total energy. From 
this point of view it is also clear that what happens at early times cannot 
possibly affect the late time radiation, since the produced pairs only see the 
late time geometry. The conventional derivation of radiance obscures this 
point somewhat and it seems desirable to find an approach which makes 
this feature manifest from the outset. For the two systems considered here, 
and presumably this is true in general, the effect of the tunneling was to 
shift the distribution of any particles that were present before tunneling. In 
the present case initial excitations were damped because the final surface is 
rotated clockwise relative to the initial surface. A counterclockwise rotation 
would have led to amplification. In |48| numerical investigations are quoted 



which show that the rotation is always clockwise for the false vacuum bubble. 
One is led to speculate whether this is a general phenomenon — whether all 
tunneling transitions lead to damping. 
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